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This topical review article gives an overview of the interplay between quantum information the¬ 
ory and thermodynamics of quantum systems. We focus on several trending topics including the 
foundations of statistical mechanics, resource theories, entanglement in thermodynamic settings, 
fluctuation theorems and thermal machines. This is not a comprehensive review of the diverse held 
of quantum thermodynamics; rather, it is a convenient entry point for the thermo-curious infor¬ 
mation theorist. Furthermore this review should facilitate the unification and understanding of 
different interdisciplinary approaches emerging in research groups around the world. 
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I. INTRODUCTION 

If physical theories were people, thermodynamics 
would be the village witch. Over the course of three 
centuries, she smiled quietly as other theories rose and 
withered, surviving major revolutions in physics, like the 
advent of general relativity and quantum mechanics. The 
other theories find her somewhat odd, somehow different 
in nature from the rest, yet everyone comes to her for 
advice, and no-one dares to contradict her. Einstein, for 
instance, called her “the only physical theory of universal 
content, which I am convinced, that within the frame¬ 
work of applicability of its basic concepts will never be 
overthrown.” 

Her power and resilience lay mostly on her frank inten¬ 
tions: thermodynamics has never claimed to be a means 
to understand the mysteries of the natural world, but 
rather a path towards efficient exploitation of said world. 
She tells us how to make the most of some resources, like 
a hot gas or a magnetized metal, to achieve specific goals, 
be them moving a train or formatting a hard drive. Her 
universality comes from the fact that she does not try to 
understand the microscopic details of particular systems. 
Instead, she only cares to identify which operations are 
easy and hard to implement in those systems, and which 
resources are freely available to an experimenter, in order 
to quantify the cost of state transformations. Although it 
may stand out within physics, this operational approach 
can be found in branches of computer science, economics 
and mathematics, and it plays a central role in quantum 
information theory—which is arguably why quantum in¬ 
formation, a toddler among physical theories, is bringing 
so much to thermodynamics. 

In the early twentieth century, information theory was 
constructed as the epitome of detachment from physics 
[7]. Its basic premise was that we could think of infor¬ 
mation independently of its physical support: a message 
in a bottle, a bit string and a sensitive phone call could 
all be treated in the same way. This level of abstraction 
was not originally conceived for its elegance; rather, it 
emerged as the natural way to address very earthly ques¬ 
tions, such as “can I reliably send a message through a 
noisy line?” and “how much space do I need to store a 
picture?”. In trying to quantify the resources required 
by those tasks (for example, the number of uses of the 
noisy channel, or of memory bits), it soon became clear 
that the relevant quantities were variations of what is 
now generally known as entropy [8]. Entropy measures 
quantify our uncertainty about events: they can tell us 
how likely we are to guess the outcome of a coin toss, 
or the content of a message, given some side knowledge 
we might have. As such, they depend only on proba¬ 
bility distributions over those events, and not on their 
actual content (when computing the odds, is does not 
matter whether they apply to a coin toss or to a horse 
race). Information theory has been greatly successful in 
this approach, and is used in fields from file compression 
to practical cryptography and channel coding [.' ]. 


But as it turned out, not all information was created 
equal. If we zoom in and try to encode information in the 
tiniest support possible, say the spin of an electron, we 
face some of the perplexing aspects of quantum physics: 
we can write in any real number, but it is only possible 
to read one bit out, we cannot copy information, and we 
find correlations that cannot be explained by local theo¬ 
ries. In short, we could not simply apply the old informa¬ 
tion theory to tasks involving quantum particles, and the 
scattered study of quirky quantum effects soon evolved 
into the fully-fledged discipline of quantum information 
theory [9] . Today we see quantum theory as a generaliza¬ 
tion of classical probability theory, with density matrices 
replacing probability distributions, measurements taking 
the place of events, and quantum entropy measures to 
characterize operational tasks [10]. 

While quantum information theory has helped us un¬ 
derstand the nature of the quantum world, its practical 
applications are not as well spread as for its classical 
counterpart. Technology is simply not there yet—not at 
the point where we may craft, transport and preserve all 
the quantum states necessary in a large scale. These tech¬ 
nical limitations, together with a desire to pin down ex¬ 
actly what makes quantum special, gave rise to resource 
theories within quantum information, for instance the¬ 
ories of entanglement [13]. There, the rough premise is 
that entangled states are useful for many interesting tasks 
(like secret key sharing), but distributing entanglement 
over two or more agents by transporting quantum parti¬ 
cles over a distance is hard, as there are always losses in 
the process [ 4]. Therefore, all entangled states become 
a precious resource, and we study how to distill entangle¬ 
ment from them using only a set of allowed operations, 
which are deemed to be easier to implement—most no¬ 
toriously, local operations and classical communication 
[ ]■ 

Other resource theories started to emerge within 
quantum information -purity and asymmetry have also 
been framed as resources under different sets of con¬ 
straints—and this way of thinking quickly spread among 
the quantum information community (see Ref. [16] for a 
review). As many of its members have a background in 
physics and an appetite for abstraction, it was a natural 
step for them to approach thermodynamics with such a 
framework in mind. Their results strengthen thermody¬ 
namics, not only by extending her range of applicabil¬ 
ity to small quantum systems, but also by revisiting her 
fundamental principles. The resource theory approach to 
thermodynamics is reviewed in Section III. 

Each resource theory explores the limitations imposed 
by one specific physical constraint, like locality or energy 
conservation. In a realistic setting we could be bound to 
several of these constraints, a natural case that can be 
modelled by combining different resource theories, thus 
restricting the set of allowed operations. In Section IV 
we review and discuss attempts to combine thermody¬ 
namic and locality constraints. In particular, we look 
at the role of entanglement resources in thermodynamic 
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An early puzzle in thermodynamics: imagine a box filled with a gas, with a partition in the middle. An 
agent (the demon) who can observe the microscopic details of the gas particles, controls a small gate in the 
partition, selectively opening it to let slow particles flow to the left and fast ones to the right. This creates 
a temperature differential between the two sides. The demon can exploit this difference to extract work, by 
letting the hot gas on the right expand. This apparent contradiction with the second law of thermodynamics 
can be easily explained from an information-theoretical viewpoint. The demon had access to much more 
information than the standard observer assumed in the derivation of traditional thermodynamics, who can 
only read a few macroscopic parameters of a gas and assumes a uniform distribution over all compatible 
micro-states. Therefore, it seems natural that the demon may extract more work than predicted by standard 
thermodynamics—and this insight motivates the need for thinking of thermodynamics as a subjective resource 
theory, and extending it to the quantum regime. In the larger picture, Bennett showed that the amount of 
work needed to erase the demon’s memory at the end of the procedure (or equivalently, to prepare the memory 
to store the necessary information on the particles in the beginning) precisely makes up for the work extracted 
[1]. For reviews, see [2-6]. 


tasks, thermodynamic witnesses of non-classicality, and 
entanglement witnesses in phase transitions. 

Information theory also shed light on fundamental is¬ 
sues in statistical mechanics - the mathematical back¬ 
bone of thermodynamics. Perhaps one of the earliest sig¬ 
nificant contributions is the maximal entropy principle 
introduced by Jaynes [17, 18]. In these seminal works 
Jaynes addresses the issue of justifying the methods of 
statistical mechanics from microscopic mechanical laws 
(classical or quantum) using tools from information the¬ 
ory. In fact, deriving statistical mechanics, and hence 
thermodynamics from quantum mechanics is almost as 
old as quantum mechanics itself starting with the work 
of von Neumann [11, 19]. This is still very much and 
ongoing and active research area and in recent years has 
received significant attention from the quantum informa¬ 
tion community. The most significant contributions are 
reviewed in Section II. 

In the past twenty years, the field of non-equilibrium 
statistical mechanics has seen a rapid development in 
the treatments of driven classical and quantum systems 
beyond the linear response regime. This has culmi¬ 
nated in the discovery of various fluctuation theorems 
which relate equilibrium thermodynamic quantities to 


non-equilibrium ones, and led to a revision on how we un¬ 
derstand the thermodynamics of systems far from equi¬ 
librium [20-2 ]. Although this approach is relatively 
recent from a statistical physics perspective, a cross¬ 
fertilisation with concepts ubiquitous in quantum infor¬ 
mation theory has already started, including phase esti¬ 
mation techniques for extraction of work and heat statis¬ 
tics and feedback fluctuation theorems for Maxwell’s 
demons. In Section V we identify these existing rela¬ 
tionships and review areas where more overlap could be 
developed. 

As ideas and concepts emerge and develop it is not sur¬ 
prising that quantum information theorists have started 
to turn towards the pragmatic goal of describing the ad¬ 
vantages and disadvantages of machines which operate 
at and below the quantum threshold. Although ideas re¬ 
lating quantum engines have been around now for a sur¬ 
prising long time [25-27] - questions pertaining to the in¬ 
trinsic quantumness in the functioning of such machines 
have been raised using the tools of quantum information 
theory only relatively recently. We review progress along 
these lines in Section VI. In summary, we will review 
landmark and recent articles in quantum thermodynam¬ 
ics, discuss different approaches and models, and peek 















4 


FIG. 2: The thermodynamic origin of the von Neumann entropy 
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In 1932, von Neumann designed this thought experiment to determine the entropy of a density operator p [11]. 
The experiment accounts for the work cost of erasing the state of a gas of n atoms, initially in an ensemble 
p®", with p = J2 k Pk \(j>k) (</>fc|, by transforming it into a pure state |</>i) <2m by means of a reversible process. 

It consists of 3 steps: 1. Separation of the species: the atoms in different states \4>i),.. . ,\4> m ) inside a box 
of volume V are separated in different boxes of the same volume V by means of semi-permeable walls (from 
a to b and finally c). Note that no work has been done and no heat has been exchanged. 2. Compression: 
every gas \tf>k) is isothermally compressed to a volume 14 = PkV (from c to d). The mechanical work done 
in that process is Wk = npk ln(I4/T) = p^lnpk- The total entropy increase per particle of that process is 
AS = YlkP k ^ n Pk- 3. Unitary transformation: every gas is put in the | <f>\) state by applying different unitary 
transformations \(f>k) —> |<(>i), which are taken for free (from d to e). As the entropy of the final state is zero, 
the entropy of the initial ensemble reads S(p) = — Tr(plnp). 

Historically, it is remarkable that the Shannon entropy, which can be seen as a particular case of the von 
Neumann entropy for classical ensembles, was not introduced until 1948 [7], and Landauer’s principle was 
proposed only in 1961 [12]. 


into future directions of the field. 


SCOPE AND OTHER REVIEWS 

This review focuses on landmark and recent articles in 
the field of quantum thermodynamics with a special em¬ 
phasis on contributions from quantum information the¬ 
ory. We place emphasis on current trending topics, dis¬ 
cuss different approaches and models and peek into the 
future directions of the field. As the review is “topi¬ 
cal” , we intend to focus our attention on the interplay of 
quantum information and thermodynamics and we have 
written such that interested readers from different com¬ 
munities will get a detailed overview of how their respec¬ 
tive techniques have been successfully applied to provide 
a deeper understanding of the field. 

As the vastness of possible topics could easily exceed 
the scope of a topical review, we refer to other review 
articles and books concerning questions that have already 
been covered by other authors: 


• Equilibration and thermalization. Recovering 
statistical mechanics from the unitary evolution of 
a closed quantum system is an issue which is almost 
as old as quantum mechanics itself. This topic, far 
from being an academic issue, has seen an unprece¬ 
dented revival of interest due mainly to advances in 
experimental ultra-cold atoms. We discus the topic 
in Section II, from a quantum information perspec¬ 
tive. This topic is more extensively reviewed in 
Ref. [28]. For readers interested in this topic from 
a condensed matter perspective we recommend the 
review [29] and the special issue [30] for more recent 
developments. 

• Thermal machines. As mentioned in the intro¬ 
duction viewing engine cycles from a fully quantum 
mechanical perspective is also not a new topic [25- 
!7]. Many results on quantum engines exist which 
are not directly related to quantum information 
processing we exclude them from Section VI and 
the interested reader may learn more in Refs. [31- 
33]. 
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• Maxwell’s demon and Landauer’s principle. 

Almost as old as thermodynamics itself is the 
Maxwell’s demon paradox, briefly introduced in 
Figure 1 and Example 3. The demon paradox in¬ 
spired the seminal work of Szilard to reformulate 
the demon as a binary decision problem [34], The 
resolution of Maxwell demon paradox by Landauer 
cements the relationship between the physical and 
information theoretical worlds. This demon has 
been extensively investigated from both a quantum 
and classical perspective in Refs. [2-6]. 

• Quantum thermodynamics. The 2009 book 
[35] covers a range of topics regarding the emer¬ 
gence of thermodynamic behaviour in composite 
quantum systems. 

• Entanglement and phase transitions in con¬ 
densed matter. Entanglement is frequently used 
as an indicator of quantum phase transitions in con¬ 
densed matter systems. We do not cover this par¬ 
ticular setting but the interested reader may find a 
comprehensive review in Ref. [36]. 

• Resource theories. Examples and common fea¬ 
tures of resource theories (beyond quantum infor¬ 
mation theory) are discussed in Ref. [37]. In par¬ 
ticular, different approaches to general frameworks 
are discussed in Section 10 of that work. 

• Experimental implementations. Experiments 
with demons, thermal engines and work extraction 
are discussed in more depth in the perspective ar¬ 
ticle [38]. 


Definitions and notation 


Conventions followed unless otherwise stated: 
States. Discrete Hilbert spaces CU. States p are 
represented by Hermitian matrices (Tr(p) = 1 and 
p > 0). Subsystems are denoted by Roman sub¬ 
scripts, p A := Tr b {pab)- 

Entropy. Von Neumann entropy with base 2 log¬ 
arithm, S{p) = - Tr (p log 2 (p)). 

Mutual information. Measures correlations, 
I(A : B) p := S(p A ) + S(p B ) - S(p AB ). 

Energy. Hamiltonian H , average energy (H) p = 
Tr (pH), eigenvalues { E k } k , eigenvectors {| E k )} k , 
or {\E l k )}k,i if there are degeneracies, with energy 
projectors B k = | E{) (. E{\ . 

Thermal states. Gibbs state r(/3) = e „ , with 
partition function Z = Tr(e~^ H ) and inverse tem¬ 
perature /3 := -jApr- 

Free energy. F p (p) : = (H) p - ^pS(p). 
Linbladian. C(p) generates Markovian, time- 
homogeneous, non-unitary dynamics. 


II. FOUNDATIONS OF STATISTICAL 
MECHANICS 

At first sight, thermodynamics and quantum theory 
are incompatible. While thermodynamics and statisti¬ 
cal mechanics state that the entropy of the universe as a 
whole is a monotonically increasing quantity, according 
to quantum theory the entropy of the universe is constant 
since it evolves unitarily. This leads us to the question 
of to which extent the methods of statistical physics can 
be justified from the microscopic theory of quantum me¬ 
chanics and both theories can be made compatible. Un¬ 
like classical mechanics, quantum mechanics has a way 
to circumvent this paradox: entanglement. We observe 
entropy to grow in physical systems because they are en¬ 
tangled with the rest of the universe. In this section we 
review the progress made on this topic in recent years 
which show that equilibration and thermalization are in¬ 
trinsic to quantum theory. 

A. Equal a priori probabilities postulate as a 
consequence of typicality in Hilbert spaces 

Let us consider a closed system that evolves in time re¬ 
stricted to some global constraint. The principle of equal 
a priori probabilities states that, at equilibrium, the sys¬ 
tem is equally likely to be found in any of its accessible 
states. This assumption lies at the heart of statistical 
mechanics since it allows for the computation of equilib¬ 
rium expectation values by performing averages on the 
phase space. However, there is no reason in the laws 
of mechanics (and quantum mechanics) to suggest that 
the system explores its set of accessible states uniformly. 
Therefore, the equal a priori probabilities principle has 
to be put in by hand. 

One of the main insights from the field of quantum 
information theory to statistical mechanics is the substi¬ 
tution of the Equal a priori probabilities postulate by the 
use of typicality arguments [39, 40]. To be more precise, 
let us consider a quantum system described by a Hilbert 
space Hs®Hb where Hs contains the degrees of freedom 
that are experimentally accessible and Hb the ones that 
are not. In practice, we think of S as a subsystem that we 
can access, and B as its environment (sometimes called 
the bath). Concerning the global constraint, in classical 
mechanics, it is defined by the constants of motion of the 
system. In quantum mechanics, we model the restriction 
as a subspace Hr C Hs ®Hb- Let us denote by cIr , ds 
and ds the dimensions of the Hilbert spaces Hr, Hs and 
Hb respectively. 

The equal a priori probability principle would describe 
the equilibrium state as 


and would imply the state of the subsystem S to be 


Us = Trs Sr . 


( 2 ) 
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In Ref. [39] it is shown that, if we look only at the sub¬ 
system S, most of the states in Hr are indistinguishable 
from the equal a priori probability state, i. e. for most 
\ip) € Hr, Trg | ip) {4>\ ~ fls- More explicitly, if |i/>) is 
randomly chosen in Hr according to the uniform distri¬ 
bution given by the Haar measure, then the probability 
that Tr s \ip) (i/j\ can be distinguished from fig decreases 
exponentially with the dimension of Hr, cIr 

Prob[||Tr s (|V’} (ip\) - Hslli > d^ 1/3 ] <2exp (-C<4 /3 ) , 

( 3 ) 

where C is a constant and [| -1| j is the trace norm. The 
trace norm ||p — ct|| 1 measures the physical distinguisha- 
bility between the states p and a in the sense that a 
||p — o'IIi = su Po<i I Tr(Op) — Tr(Ocr)|, where the maxi¬ 
mization is made over all the observables O with operator 
norm bounded by 1. The proof of Eq. (3) relies upon con¬ 
centration of measure and in particular on Levy’s Lemma 
(see Ref. [39] for details). Let us mention that ideas in 
this spirit can be already found in S. Lloyd’s Ph.D. thesis 
[41] published in 1991. In particular, he presents bounds 
on how the expectation values of a fixed operator taken 
over random pure states of a restricted subspace fluctu¬ 
ate. 

The weakness of the previous result lies in that the use 
of typicality is made in the whole subspace Hr and, as 
we will justify next, this is not a physical assumption. In 
nature, Hamiltonians have local interactions and systems 
evolve for times that are much smaller than the age of 
the universe. Most states in the Hilbert space simply 
cannot be generated by evolving an initial product state 
under an arbitrary time-dependent local Hamiltonian in 
a time that scales polynomially in the system size [42]. 
Therefore, sampling uniformly from the whole Hilbert 
space is not physically meaningful. There has been a 
strong effort to generalize the concept of typicality for 
different sets of states [43-45]. 

The first “realistic” set of states in which typicality 
was studied was the set of matrix product states (MPS) 
[46, 47]. These type of states have been proven to de¬ 
scribe ground states of one-dimensional gapped Hamil¬ 
tonians. They are characterized by the rank of a bipar¬ 
tition of the state. This parameter quantifies the maxi¬ 
mum entanglement between partitions of an MPS. The 
MPSs with fixed rank form a set of states with an efficient 
classical representation (they only require polynomial re¬ 
sources in the number of particles). In Ref. [43], it is 
proven that typicality occurs for the expectation value of 
subsystems observables when the rank of the MPS scales 
polynomially with the size of the system with a power 
greater than 2. 

Another set recently considered in the literature has 
been the so called set of physical states which consists 
of all states that can be produced by evolving an ini¬ 
tial product state with a local Hamiltonian for a time 
polynomial in the number of particles n. By Trotter de¬ 
composing the Hamiltonian, such a set can be proven to 


be equivalent to the set of local random quantum circuits, 
that is, quantum circuits of qubits composed of polyno¬ 
mially many nearest neighbour two-qubit gates [42]. In 
Ref. [48], it was shown that the local random quantum 
circuits form an approximate unitary 2-design, i. e. that 
random circuits of only polynomial length will approx¬ 
imate the first and second moments of the Haar distri¬ 
bution. In Ref. [44] the previous work was extended to 
poly(n)-designs. Finally, let us mention that the entan¬ 
glement properties of typical physical states were studied 
in Ref. [15]. 

Let us mention that /c-designs also appear naturally in 
the context of decoupling theorems in which a the sub¬ 
system S undergoes a physical evolution separated from 
the environment B, and one wonders under what con¬ 
ditions this evolution destroys all initial correlations be¬ 
tween S and B. In particular, in Ref. [49] it is shown 
that almost-2-designs decouple the subsystem S from B 
independently of B 's size. 

Another objection against typicality is that there are 
many physically interesting systems, e. g. integrable mod¬ 
els, which, although their initial state belongs to a cer¬ 
tain restricted subspace Hr, their expectation values dif¬ 
fer from the completely mixed state in R, £r, as ex¬ 
pected from typicality arguments. This is a consequence 
of the fact that their trajectories in the Hilbert subspace 
Hr don’t lie for the overwhelming majority of times on 
generic states (see Fig. 3). Hence, in practice, state¬ 
ments on equilibration and tlrermalization will depend 
on the dynamical properties of every system, that is, on 
their Hamiltonian. This leads us to the notion of dy¬ 
namical typicality. In contrast to the kinematic typical¬ 
ity presented in this section, where an ensemble has been 
defined by all the states that belong to a certain sub¬ 
space, in dynamical typicality the ensemble is defined by 
all states that share the same constants of motion given 
a Hamiltonian H and an initial state |^>(0)). Studying 
whether typicality also holds in such a set will be pre¬ 
cisely the problem addressed in the next section. 


B. Equilibration. Maximum entropy principle from 
quantum dynamics 

In this context of deriving thermodynamics from quan¬ 
tum mechanics the first problem that needs to be ad¬ 
dressed is equilibration, that is, understand how the re¬ 
versible unitary dynamics of quantum mechanics make 
systems equilibrate and evolve towards a certain state 
where they remain thereafter. 

Because of the unitary dynamics, equilibration is only 
possible if the set of observables is restricted. In this 
spirit, a set of sufficient conditions for equilibration to¬ 
wards the time averaged state has been presented for lo¬ 
cal observables [50, 5 ] and observables of finite precision 
[52, 53]. The two approaches are proven to be equivalent 
in Ref. [54] and it is remarkable that the conditions given 
are weak and naturally fulfilled in realistic situations. 
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Scheme of the restricted subspace Hr with its untypical states forming little islands coloured in yellow. The 
left trajectory (dashed line) passes mostly on typical states while the right trajectory (solid line) has a non- 
negligible support on states that are not typical. 


For simplicity, let us here focus on equilibration of 
subsystems and, as above, identify in the total system 
a subsystem S and its environment B. The dynam¬ 
ics of the total system are governed by the Hamilto¬ 
nian H with eigenvalues {Ek}k and eigenvectors {\Ek)}k- 
This leads to the time evolution \4>(t)) = e~ lHt |^(0)) 
and the reduced state of S is psit) = Ti 'b p{t) with 
p{t) = \ip(t))(ip(t)\. 

If equilibration happens, then it happens towards the 
time averaged state i. e. u>s '■= Tr^w with 

u = lim ^ [ p(t)dt = Y Pk P (0)Pk (4) 
T^oo T J 0 ^ 

with Pk the projectors onto the Hamiltonian eigenspaces. 
The time averaged state is the initial state dephased in 
the Hamiltonian eigenbasis. For this reason it is also 
called diagonal ensemble. 

In Ref. [50], a notion of equilibration is introduced by 
means of the average distance (in time) of the subsys¬ 
tem psit) from equilibrium. A subsystem S is said to 
equilibrate if 

(\\Ps{t) - 11x>t := lim ^ / dt ||p s (t) - < 1, 

T-^oo 1 J 0 

( 5 ) 

where ||ps(t) — w s|li is the trace distance. If this average 
trace distance can be proven to be small, then the sub¬ 
system S is indistinguishable from being at equilibrium 
for almost all times. 

Equilibration as a genuine property of quantum me¬ 
chanics is shown in Ref. [50] by precisely proving that this 
average distance is typically small. More concretely, if 
the Hamiltonian that dictates the evolution of the system 
has non-degenerate gaps i. e. all the gaps of the Hamilto¬ 
nian are different (an assumption which we will comment 


on below), then the average distance from equilibrium is 
bounded by 

<IM*) - wslli)* < ^ ]id""\u;) ’ (6) 

where d eS {p) := 1/Tr(p 2 ) is the effective dimension of 
p and lob = Trsw. Roughly speaking, the effective di¬ 
mension of a state tells us how many eigenstates of the 
Hamiltonian support such state. It can also be related to 
the 2-Renyi entanglement entropy by S 2 ip) = log d eff (p). 
Hence, equation (6) guarantees equilibration for Hamil¬ 
tonians with non-degenerate energy gaps as long as the 
initial state is spread over many different energies. 

Although the condition of having non-degenerate gaps 
may look very restrictive at first sight, note that Hamilto¬ 
nians that do not fulfil it form a set of zero measure in the 
set of Hamiltonians, since any arbitrarily weak perturba¬ 
tion breaks the degeneracy of the gaps. In Ref. [51], the 
non-degenerate gaps condition was weakened by showing 
that equilibration occurs provided that no energy gap is 
hugely degenerate. This condition can be understood as 
a way of preventing the situation where there is a sub¬ 
system which does not interact with the rest. 

Let us finally point out that the equilibrium state intro¬ 
duced in Eq. (4) is precisely the state that maximizes the 
von Neumann entropy given all the conserved quantities 
[55]. This observation turns the principle of maximum 
entropy into a consequence of the quantum dynamics. 
The principle of maximum entropy was introduced by 
Jaynes in Ref. [17] and states that the probability distri¬ 
bution which best represents the current state of knowl¬ 
edge of the system is the one with largest entropy given 
the conserved quantities of the system. We will come 
back in more detail to the Jaynes principle in the next 
section when the thermalization for integrable systems is 












discussed. 


C. Thermalization. Emergence of Gibbs states in 
local Hamiltonians 

The next step in this program of justifying the meth¬ 
ods of statistical mechanics from quantum mechanics is 
to tackle the issue of thermalization, i. e. to understand 
why the equilibrium state is usually well described by a 
Gibbs state, which is totally independent of the initial 
state of the system, except for some macroscopic con¬ 
straints such as its mean energy. In Ref. [56], a set of 
sufficient conditions for the emergence of Gibbs states 
is presented for the case of a subsystem S that inter¬ 
acts weakly with its environment B through a coupling 
V. The Hamiltonian that describes such a situation is 
H = Hs + Hb + V. These conditions are a natural 
translation of the three ingredients that enter the stan¬ 
dard textbook proof of the canonical ensemble in classical 
statistical physics: 

1. The equal a priory probability postulate that has 
been replaced by typicality arguments in Sec¬ 
tion II A, and an equilibration postulate (such as 
the second law) that has been replaced by quan¬ 
tum dynamics in Section IIB. 

2. The assumption of weak-coupling. Here, the stan¬ 
dard condition from perturbation theory, HRH^ -C 
gaps {H), is not sufficient in the thermodynamic 
limit, due to the fast growth of the density of states 
and the corresponding shrinking of the gaps in the 
system size. Instead, it is replaced with a physi¬ 
cally relevant condition, HVjl^ <C fee T, which is 
robust in the thermodynamic limit. 

3. An assumption about the density of states of the 
bath [57] , namely, that it grows faster than expo¬ 
nentially with the energy and that it can be locally 
approximated by an exponential. 

Note that the weak-coupling condition will not be sat¬ 
isfied in spatial dimensions higher than one for suffi¬ 
ciently large subsystems, since the interaction strength 
typically scales as the boundary of the subsystem S. This 
will be the case regardless of the strength of the coupling 
per particle or the relative size between S and B. This 
should not be seen as a deficiency of the above results, 
but as a feature of strong interactions. Systems that 
strongly interact with their environment do not in gen¬ 
eral equilibrate towards a Gibbs state, in a similar way 
that the reduced state (of a part) of a Gibbs state need 
not have Gibbs form [58, 59]. In this context, the findings 
of Ref. [60] suggest that subsystems do not relax towards 
a local Gibbs state but to the reduction of the global 
Gibbs state; this is shown for translation-invariant quan¬ 
tum lattices with finite range but arbitrarily strong inter¬ 
actions. The eigenstate thermalization hypothesis (ETH) 


[61, 62] gives further substance to this expectation. ETH 
has several formulations. Its simplest one is maybe the 
one introduced in Ref. [62]. It states that the expectation 
value (Ek\0\Ek) of a few-body observable O in an in¬ 
dividual Hamiltonian eigenstate | Ek) equals the thermal 
average of O at the mean energy Ek . Although ETH has 
been observed for some models, it is not true in general 
and it is well known to break down for integrable mod¬ 
els (see [62] for an example with hard-core bosons and 
references in Ref. [28] for further examples). 

In the same spirit, it has recently been proven that a 
global microcanonical state (the completely mixed state 
of a energy shell subspace spanned by the Hamiltonian 
eigenstates with energy inside a narrow interval) and a 
global Gibbs state are locally indistinguishable for short 
range spin Hamiltonians off criticality, that is, when they 
have a finite correlation length [63]. This represents a 
rigorous proof of the so called equivalence of ensembles. 
If the Hamiltonian is not translationally invariant, the 
local indistinguishability between canonical and micro- 
canonical ensembles becomes a typical property of the 
subsystems, allowing for rare counterexamples. 

Concerning the latter condition on the density of states 
of the bath, in Ref. [64] it has been proven that the den¬ 
sity of states of translational invariant spin chains tends 
to a Gaussian in the thermodynamic limit, matching the 
suited property of being well approximated by an expo¬ 
nential. In Ref. [63], the same statement is proven for 
any short ranged spin Hamiltonian. 

Let us finally point out that not all systems thermalize. 
For instance, integrable systems are not well described by 
the Gibbs ensemble. This is due to the existence of local 
integrals of motion, i. e. conserved quantities, Q a that 
keep the memory about the initial state. Instead, they 
turn to be described by the Generalized Gibbs Ensemble 
(GGE) defined as 

TGGE oc exp ^-/3 ^ H + ^ ^ (7) 

where the generalized chemical potential p a is a Lagrange 
multiplier associated to the conserved quantity Q a such 
that its expectation value is the same as the one of the ini¬ 
tial state. The GGE was introduced by Jaynes in Ref. [17] 
where he pointed out that statistical physics can be seen 
as statistical inference and an ensemble as the least bi¬ 
ased estimate possible on the given information. Nev¬ 
ertheless, note that any system has as many conserved 
quantities as the dimension of the Hilbert space, e.g. 
Qa = | E a ) (E a \. If one includes all these conserved quan¬ 
tities into the GGE the ensemble obtained is the diagonal 
ensemble introduced in Eq. (4). Note that the descrip¬ 
tion of the equilibrium state by the diagonal ensemble 
requires the specification of as many conserved quanti¬ 
ties as the dimension of the Hilbert space, which scales 
exponentially in the system size, and becomes highly in¬ 
efficient. A question arises here naturally, is it possible to 
provide an accurate description of the equilibrium state 
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specifying only a polynomial number of conserved quan¬ 
tities? If so, what are these relevant conserved quantities 
Q a that allow for an accurate and efficient representation 
of the ensemble? This question is tackled in Ref. [65]. 
There, it is argued that the relevant conserved quanti¬ 
ties are the ones that make the GGE as close as possible 
to the diagonal ensemble in the relative entropy distance 
D(ui\\tgge), which in this particular case can be written 
as 

D(uj\\tgge) = S(tgge) ~ S(u ), (8) 

where we have used that the diagonal ensemble and 
the GGE have by construction the same expectation 
values for the set of selected conserved quantities, 
i. e. Tr(Q a TGGE) = T i(Q a uj). Equation (8) tells us that 
the relevant conserved quantities are the ones the min¬ 
imize the entropy S(tgge )• Note that in contrast to 
Jaynes approach, where entropy is maximized for a set of 
observables defined beforehand, here the notion of phys¬ 
ically relevant is provided by how much an observable is 
able to reduce the entropy by being added into the set of 
observables that defines the GGE. 

If instead of calculating the relative entropy between 
the diagonal ensemble and the GGE’s we do it with re¬ 
spect to the set of product states, i. e. 

T(u>) := min D(u>\\iri (g> 7 t 2 (g> ... 7r n ), (9) 

7T i ,7T2 , • •. ,7V n 

then we obtain a measure of the total (multipartite) corre¬ 
lations of the diagonal ensemble. In Ref. [66] the scaling 
with system size of the total correlations of the diagonal 
ensemble has been shown to be connected to ergodic- 
ity breaking and used to investigate the phenomenon of 
many-body localization. 


D. Equilibration times 

Maybe the major challenge that is still open in the 
equilibration problem is to determine the equilibration 
timescale. It turns out that even if we know that a system 
equilibrates, there are no relevant bounds on how long 
the equilibration process takes. There could be quan¬ 
tum systems that are going to equilibrate, but whose 
equilibration times are of the order of magnitude of the 
age of the universe, or alternatively, some systems, like 
glasses, which do not relax to equilibrium at all, but have 
metastable states with long lifetimes. The problem of es¬ 
timating equilibration timescales is thus essential in order 
to have a full understanding of thermalization. 

So far, progress on this issue has taken place from two 
different approaches. On the one hand, rigorous and 
completely general bounds on equilibration times have 
been presented in Ref. [51]. Due to their generality, 
these bounds scale exponentially with the system size, 
leading to equilibration times of the age of the universe 
for macroscopic systems. On the other hand, very short 


equilibration times have been proven for generic observ¬ 
ables [67], Hamiltonians [68-72], and initial states [73]. 
In nature, systems seem to equilibrate in a time that is 
neither microscopic nor exponential in the system size. A 
relevant open question is what properties of the Hamilto¬ 
nians and operators lead to reasonable equilibration time. 
As a first step, in Ref. [74], a link between the complexity 
of the Hamiltonian’s eigenvectors and equilibration time 
is presented. The result does not completely solve the 
question, since the given bounds are not fulfilled by all 
Hamiltonians but only by a fraction of them, and further 
research in this direction is needed. 


E. Outlook 

The aim of this section has been to justify that thermal 
states emerge in Nature for generic Hamiltonians. To 
complete the picture presented here we recommend the 
article [28] where an extensive review of the literature on 
foundations of statistical mechanics is provided. 

The main ideas presented here have also been widely 
studied in the context of condensed matter physics, in 
which systems are typically brought out of equilibrium 
by sudden (and slow) quantum quenches: the Hamilto¬ 
nian of a system (that is initially in the ground state) 
is suddenly (or smoothly) changed in time. We recom¬ 
mend the review article [29] on non-equilibrium dynamics 
of closed interacting quantum systems. 

Let us finish the section with a list of some of the open 
problems that we consider most relevant in the field: 

• Typicality for symmetric states. Hamiltonians in 
nature are not generic but have symmetries. Hence, 
the notion of typicality should be extended to phys¬ 
ical states that are produced by symmetric Hamil¬ 
tonians. 

• Quantum notion of integrability. One of the reasons 
why it is so difficult to extract strong statements 
on the equilibration and thermalization of many 
body quantum systems is the absence of a satis¬ 
factory quantum notion of integrability [75]. This 
leads first to some widespread confusion, since in¬ 
tegrability is mentioned very often in the field of 
non-equilibrium dynamics, and second it does not 
allow us to classify quantum systems into classes 
with drastically different physical behaviour, like 
what occurs in classical mechanics. 

• Equilibration times. Without bounds on the equi¬ 
libration time scales, statements on equilibration 
become useless. As we have seen, the equilibration 
times are model dependent. We need then to un¬ 
derstand how the equilibration times depend on the 
features of the Hamiltonian and the set of observ¬ 
ables considered. 

• Relative thermalization. It was highlighted in 
Ref. [76] that local thermalization of a subsystem S 7 
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as described here, is not enough to guarantee that 
S will act as thermal bath towards another physical 
system R. In other words, imagine that we want to 
perform quantum thermodynamics on a reference 
system R , using S as a thermal bath. To model a 
thermodynamic resource theory that recovers the 
laws of thermodynamics, it is not sufficient to de¬ 
mand that S be in a local Gibbs state rg(/3). Indeed 
we need S to be thermalized relative to R , that 
is the the two systems should be uncorrelated, in 
global state, Ts(f3)® Pr. If this does not hold, then 
we cannot recover the usual thermodynamic mono- 
tones (for instance, there could be anomalous heat 
flows against the temperature gradient). Therefore, 
the relevant question for resource theories of ther¬ 
modynamics is not only “does S thermalize locally 
after evolving together with an environment?”, but 
rather ‘does S thermalize relative to R after evolv¬ 
ing together with an environment?’, and the results 
discussed in this section should be generalized to 
that setting. First steps in this direction can be 
found in Ref. [76], where the authors use decou¬ 
pling —a tool developed in quantum information 
theory—to find initial conditions on the entropies 
of the initial state that lead to relative thermaliza- 
tion. 


III. RESOURCE THEORIES 

In the previous section we saw the progress that has 
been made in understanding how systems come to equi¬ 
librium, in particular thermal equilibrium, and as such 
a justification for the thermal state. In the rest of this 
review we will now take the thermal state as a given, 
and see what is the thermodynamics of quantum systems 
which start thermal or interact with thermal states. We 
will start from an operational point of view, treating the 
thermal state as a “free resource”, a view inspired by 
other resource theories from quantum information. 

In this section we discuss the approach of thermo¬ 
dynamics as a resource theory in more detail. Let us 
start by introducing the basic ideas behind resource the¬ 
ories that can be found in the literature, entanglement 
theory being the paradigmatic example. The first step 
is to fix the state space S, which is usually compati¬ 
ble with a composition operation—for instance, quan¬ 
tum states together with the tensor product, in systems 
with fixed Hamiltonians. The next step is to define the 
set of allowed state transformations. For thermodynam¬ 
ics, these try to model adiabatic operations- like energy¬ 
preserving reversible operations, and contact with a heat 
bath. 

The set of allowed operations induces a structure on 
the state space: we say that p —> a if there is an allowed 
transformation from p to a. The relation —> is a pre¬ 
order , that is, a binary relation that is both reflexive 
(p —> a) and transitive (p —> a and a —> r implies p —► 


t; this results from composing operations one after the 
other). 

The task now is to find general properties of this struc¬ 
ture. A paradigmatic example is looking for simple nec¬ 
essary and sufficient conditions for state transformations. 
The most general case are functions such that 

• p —>• a => /(p, <j ) > 0 (that is, f{p,cr) > 0 is a 
necessary condition for state transformations), or 

• f(p,a) > 0 => p —> a (that is, f(p,cr) > 0 is a 
sufficient condition for state transformations). 

Often, we try to find necessary and sufficient conditions 
as functions that can be written like f(p, a) = g(p) — 
h(a). In the special case where g = h for a necessary 
condition {p —► a => g(p) > g(cr)), we call g a monotone 
of the resource theory. For example, in classical, large- 
scale thermodynamics, the free energy is a monotone. 

In order to quantify the cost of state transformations, 
we often fix a minimal unit in terms of a standard resource 
that can be composed. For example, in entanglement 
theory the standard resource could be a pair of maximally 
entangled qubits, and in quantum thermodynamics we 
could take a single qubit (with a fixed Hamiltonian) in a 
pure state. The question then is ‘how many pure qubits 
do I need to append to p in order to transform it into 
<r?’ or, more generally, ‘what is the cost or gain, in terms 
of this standard resource, of the transformation p —> oV 
[77-79], 

One may also try to identify special sets of states. The 
most immediate one would be the set of free states: those 
that are always reachable, independently of the initial 
state. In standard thermodynamics, these tend to be 
what we call equilibrium states, like Gibbs states. An¬ 
other interesting set is that of catalysts, states that can 
be repeatedly used to aid in transformations. We will 
revisit them shortly. 

A. Models for thermodynamics 

Now that we have established the basic premise and 
structure of resource theories, we may look at different 
models for resource theories of thermodynamics, which 
vary mostly on the set of allowed operations. In the good 
‘spherical cow’ tradition of physics, the trend has been to 
start from a very simple model that we can understand, 
and slowly expand it to reflect more realistic scenarios. 
In general there are two types of operations allowed: con¬ 
tact with a thermal bath and reversible operations that 
preserve some thermodynamic quantities. Each of those 
may come in different flavours. 

1. Noisy and unital operations 

In the simplest case, all Hamiltonians are fully degen¬ 
erate, so thermal states of any temperature are just fully 
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mixed states, and there are no special conserved quan¬ 
tities. In this setting, thermodynamics inherits directly 
from the theory of noisy operations [80]. We may model 
contact with a thermal bath as composition with any sys¬ 
tem in a fully mixed state, and reversible operations as 
any unitary operation. Furthermore, we assume that we 
can ignore, or trace out, any subsystem. Summing up, 
noisy operations have the form 


results are only known for classical states (states that 
are block-diagonal in the energy eigenbasis) [83] and any 
state of a single qubit [90, 91]. In the limit of a fully 
degenerate Hamiltonian, we recover the resource theory 
of noisy operations. 

3. Gibbs-preserving maps 


T(pa)= 



PA <8 


W\. 



where A! is any subsystem of AB and U is a unitary 
matrix. Alternatively, we may allow only for maps that 
preserve the fully mixed state, Ta^b '■ Ta->b(|a[) = pftj 
called unital maps (an example would be applying one 
of two isometries and then forgetting which one). The 
two sets—noisy operations and unital maps- induce the 
same pre-order structure in the state space. In this set¬ 
ting, majorization is a necessary and sufficient condition 
for state transformations [80]. Roughly speaking, ma¬ 
jorization tells us which state is the most mixed. Let 
r = (rq, r 2 ,... ,r N ) and s = (si, s 2 ,s N ) be the eigen¬ 
values of two states p and a respectively, in decreasing 
order. We say that r majorizes s if Yji=i r i > X?=i s f, 
for any k < N. In that case p —» a; monotones 
for this setting are called Schur monotone functions , of 
which information-theoretical entropy measures are ex¬ 
amples [79, 81-84]. For example, if p majorizes a, then 
the von Neumann entropy of p , S(p) = — Tr(plog 2 p), is 
smaller than S(a). For a review, see [84]. 


Following the example of the theory of noisy opera¬ 
tions, we could try to replace these thermal operations 
with so-called Gibbs-preserving maps, that is, maps such 
that Ta-+b{ta{0)) = t b(/ 3)- This constraint is easier 
to tackle mathematically, and the two resource theo¬ 
ries induce the same pre-order on classical states, lead¬ 
ing to a condition for state transformation called Gibbs- 
majorization (which is majorization after a rescaling of 
the eigenvalues) [8- ]. However, Gibbs-preserving maps 
are less restrictive than thermal operations for general 
quantum states [93]. For example, suppose that you 
have a qubit with the Hamiltonian H = E 11) (11, and 
you want to perform the transformation |1) —► |+) = 
(|0) + 11 ))/a/ 2- This is impossible through thermal op¬ 
erations, which cannot create coherence; yet there exists 
a Gibbs-preserving map that achieves the task. We may 
still use Gibbs-preserving maps to find lower bounds on 
performance, but at the moment we cannot rely on them 
for achievability results, as they are not operationally de¬ 
fined. 


4- Coherence 


2. Thermal operations 

The next step in complexity is to let systems have non¬ 
degenerate Hamiltonians. The conserved quantity is en¬ 
ergy, and equilibrium states are Gibbs states of a fixed 
temperature T. For instance for a system A with Hamil¬ 
tonian Ha, the equilibrium state is ta(P) = e~^ HA /Z. 
We can model contact with a heat bath as adding any 
system in a Gibbs state—this corresponds to the ideal¬ 
ization of letting an ancilla equilibrate for a long time. 
A first approach to model physical reversible transfor¬ 
mations is to allow for unitary operations U that pre¬ 
serve energy—either absolutely ([U,H\ = 0) or on aver¬ 
age (Tr (Hp) = Tt(H ( UpW)) for specific states). Fi¬ 
nally, we are again allowed to forget, or trace out, any 
subsystem. Together, these transformations are called 
thermal operations, 

T(pa) = T?A' (u A b [PA <8 t b (/ 3)] C/^ B j , 

where A' is any subsystem of AB and U is an energy- 
conserving unitary [85]. The monotones found so far are 
different versions of the free energy, depending on the 
exact regime [83, 86-89] (see Example 1). It is worth 
mentioning we can build necessary conditions for state 
transformations from these monotones, but sufficiency 


The difference between thermal operations and Gibbs- 
preserving maps is not the only surprise that quantum 
coherence had in store for thermodynamics enthusiasts. 
The question of how to create coherence in the first place 
led to an intriguing discovery. In order to achieve the 
above transformation |1) —> |+) through thermal oper¬ 
ations, we need to draw coherence from a reservoir. A 
simple example of a coherence reservoir would be a dou¬ 
bly infinite harmonic oscillator, H = nA |n) (n|, 

in a coherent state like |'F) = TV -1 Y^n=a \ n )• Lasers ap¬ 
proximate such reservoirs, which explains why we can use 
them to apply arbitrary transformations on quantum sys¬ 
tems like ion traps. One may ask what happens to the 
reservoir after the transformation: how much coherence 
is used up? Can we use the same reservoir to perform 
a similar operation in a new system? The unexpected 
answer is that coherence is, in a sense, catalytic: while 
the state of the reservoir is affected, its ability to imple¬ 
ment coherent operations is not [94]. What happens is 
that the state of the reservoir ‘spreads out’ a little with 
each use, but the property that determines the efficacy of 
the reservoir to implement operations stays invariant. In 
more realistic models for coherence reservoirs, where the 
Hamiltonian of the reservoir has a ground state, the cat¬ 
alytic properties hold for some iterations, until the state 
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Example 1: Free energy as a monotone. 


This is an example of finding monotones for the resource theory of thermal operations [86]. We are interested 
in finding the optimal rates of conversion between two states p and a, in the limit of many independent copies, 

R(p -»• a) := sup lim p®” -► <j® Rn . 

n—too 

If both R(p —> a),R(a —> p) > 0, and these quantities represent optimal conversion rates, then the process 
must be reversible, that is, R(p —► a) = l/i?(cr —> p); otherwise we could build a perpetual motion engine, 
and the resource theory would be trivial. The idea is to use a minimal, scalable resource a as an intermediate 
step. We can think of a as a currency: we will sell n copies of p for a number of coins, and use them to buy 
some copies of er. To formalize this idea, we define the selling and buying cost of a state p, or more precisely 
the distillation and formation rates, 

R D (p) := R{p^ a), R F (p) := R(a -S> p) = 

In the optimal limit we have the process 

p n a nR °(p) (j nRD (p) rF M ^ p n a nR (p^>°) ( 


which gives us the relation 


R(p —> er) 


R D {p) 

R D (aY 


We have reduced the question to finding the distillation rate, which depends on the choice of a. For example, 
take p, a and a to be classical states (diagonal in the energy basis) of a qubit with Hamiltonian H = A |1) (1|. 
For the currency, we choose a = |1) (1|. The distillation rate is found by use of information-compression tools 
>0]. It is given by the relative entropy between p and the thermal state r(/3), 

R D (p) =D{p\\rm 

— Tr(p(logp — logr(/3))) 

=/3(Fp(p) - F p (tU 3))), 

where Fp(p) = (E) p — /3~ 1 S{p) is the free energy of p at inverse temperature /3. All in all, we find the conversion 
rate 


R(p —► a) 


Fp{p) -Fp(r(/3)) 
Ff}(u) - Fp(T(P))' 


Now we can apply this result to find a monotone for a single-shot scenario: in order to have p —> a we need 
in particular that R(p —> a) > 1. In other words, we require Fp(p) > Fp(p), thus recovering the free energy 
as a monotone for the resource theory of thermal operations. If we work directly in the single-shot regime, we 
recover a whole family of monotones [ ] based on quantum Renyi relative entropies [92], of which the free 

energy is a member. 


spreads all the way down to the ground state. At that 
stage, the reservoir needs to be recharged with energy 
to pump up the state again. Crucially, we do not need 
to supply additional coherence. In the converse direc¬ 
tion, we know that coherence reservoirs only are critical 
in the single-shot regime of small systems. Indeed, in the 
limit of processing many copies of a state simultaneously, 
the work yields of doing it with and without access to a 
coherence reservoir converge [95] 


5. Catalysts 

The catalytic nature of coherence raises more general 
questions about catalysts in thermodynamics. Imagine 
that we want to perform a transformation p —^ er in a 
system S , and we have access to an arbitrary ancilla in 
any desired state 7 . Now suppose that our constraint is 
that we should return the ancilla in a state that is e-close 
to 7 : 

Ps®'YA-> <rsA- \Wa ~ 7a||i < e- 

The question is whether we can overcome the usual limits 
found in thermal operations by use of this catalyst. In 
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other words, can we perform the above transformation 
in cases where p —» er would not be allowed? It turns 
out that if no other restrictions are imposed on the cata¬ 
lyst, then for any finite e and any two states p and a, we 
can always find a (very large) catalyst that does the job 
[83]. These catalysts are the thermodynamic equivalent 
of embezzling states in LOCC [96]. However, if we impose 
reasonable energy and dimension restrictions on the cat¬ 
alyst, we recover familiar monotones for state transfor¬ 
mations [83, 97]. These restrictions and optimal catalysts 
result from adapting the concept of trumping relations on 
embezzling states [98, 99] to the thermodynamic setting. 
In particular, if we demand that e oc n -1 , where n is 
the number of qubits in the catalyst, we recover the free 
energy constraint for state transformations [97]. A rele¬ 
vant open question, motivated by the findings of catalytic 
coherence, is what happens if we impose operational con¬ 
straints on the final state of the catalyst. That is, instead 
of asking that it be returned e-close to 7 , according to the 
trace distance, we may instead impose that its catalytic 
properties stay unaffected. It would be interesting to see 
if we recover similar conditions for allowed transforma¬ 
tions under these constraints. 


6. Clocks 

All of resource theories mentioned allow for energy¬ 
preserving unitary operations to be applied for free. That 
is only the ‘first order’ approach towards an accurate 
theory of thermodynamics, though. Actually, in order 
to implement a unitary operation, we need to apply a 
time-dependent Hamiltonian to the systems involved. To 
control that Hamiltonian, we require very precise time¬ 
keeping—in other words, precise clocks, and we should 
account for the work cost of using such clocks. Fur¬ 
thermore, clocks are clearly out of equilibrium, and us¬ 
ing them adds a source of free energy to our systems. 
Including them explicitly in a framework for work ex¬ 
traction forces us to account for changes in their state, 
and ensures that we do not cheat by degrading a clock 
and drawing its free energy. First steps in this direc¬ 
tion can be found in [67]. There, the goal is to im¬ 
plement a unitary transformation in a system S, using 
a time-independent Hamiltonian. For this, the authors 
introduce an explicit clock system C hat runs continu¬ 
ously, as well as a weight W that acts as energy and 
coherence reservoir. The global system evolves under a 
time-independent Hamiltonian, designed such that the 
Hamiltonian applied on S depends on the position of 
the clock—which effectively measures time. The authors 
show that such a construction allows us to approximately 
implement any unitary operation on S, while still obeying 
the first and second laws of thermodynamics. Further¬ 
more, the clock and the weight are not degraded by the 
procedure (just like for catalytic coherence). In particu¬ 
lar, this result supports the idea behind the framework of 
thermal operations: that energy-conserving unitaries can 


approximately be implemented for free (if we neglect the 
informational cost of designing the global Hamiltonian). 
Note that this is still an idealized scenario, in which the 
clock is infinite-dimensional and moves like a relativistic 
particle (the Hamiltonian is proportional to the parti¬ 
cle’s momentum). A relevant open question is whether 
there exist realistic systems with the properties assigned 
to this clock, or alternatively how to adapt the protocol 
to the behaviour of known, realistic clocks. That direc¬ 
tion of research can be related to the resource theory of 
quantum reference frames [90, 100-102], An alternative 
direction would be to ask what happens if we do not have 
a clock at all—can we extract all the work from a quan¬ 
tum state if we are only allowed weak thermal contact? 
This question is studied (and answered in the negative, 
for general states) in Ref. [103]. 

7. Free states and passivity 

It is now time to question the other assumption behind 
the framework of thermal operations: that Gibbs states 
come for free. There are two main arguments to support 
it: firstly, Gibbs states occur naturally under standard 
conditions, and therefore are easy to come by; secondly, 
they are useless on their own. The first point, typicality 
of Gibbs states, is essentially the fundamental postulate 
of statistical mechanics: systems equilibrate to thermal 
states of Gibbs form. This assumption is discussed and 
ultimately justified from first principles in Section II The 
second point is more subtle. Pusz and Woronowicz first 
introduced the notion of passive states, now adapted to 
the following setting [108-110]. Let S' be a system with a 
fixed Hamiltonian H , in initial state p. We ask whether 
there is a unitary U that decreases the energy of S, that 
is 

Tr(pH) > TiiUpU^H). 

If we can find such a unitary, then we could extract work 
from S by applying U and storing the energy difference 
in a weight system. If there is no U that achieves the 
condition above, then we cannot extract energy from p, 
and we say that the state is passive. The latter applies to 
classical states (i.e., diagonal in the energy basis) whose 
eigenvalues do not increase with energy. However, sup¬ 
pose that now we allow for an arbitrary number n many 
copies of p and a global unitary U g \. The question be¬ 
comes whether 

Tr(piL) > ^ Tr(t/ gl p®” H g \), 

where H g \ is the global Hamiltonian, which is the sum of 
the independent local Hamiltonians of every system. If 
this is not possible for any n, we say that p is completely 
passive , and it turns out that only states of Gibbs form, 
p = r(/3) are completely passive. Moreover, Gibbs states 
are still completely passive if we allow each of the n sub¬ 
systems to have a different Hamiltonian, as long as all 
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Example 2: Heat engines 


The extreme case where one of our resources is in itself a second heat bath is of particular interest. This is 
a very natural scenario in traditional thermodynamics: steam engines used a furnace to heat a chamber, and 
exploit the temperature difference to the cooler environment. The study of this limit led to landmark findings 
like trains, fridges and general heat engines, and to theoretical results on the efficiency of such engines. One 
might wonder whether these findings can also be applied at the quantum scale, and especially to very small 
systems composed only of a couple of qubits [25, 104]. The answer is yes: not only is it possible to build 
two-qubit heat engines, but they achieve Carnot efficiency [105, 106]. It is possible to build heat engines that 
do not require a precise control of interactions, in other words, that do not require a clock [105, 107]. 


the states correspond to the same inverse temperature 
ft. This justifies the assumption that we may bring in 
any number and shape of subsystems in thermal states 
for free, because we could never extract work from them 
alone—another resource is necessary, precisely a state out 
of equilibrium. More formally, it was shown that if a re¬ 
source theory allows only for energy-conserving unitaries 
and composition with some choice of free states, Gibbs 
states are the only choice that does not trivialize the the¬ 
ory [83, 111]. 

8. Different baths 

The results outlined above suggest that thermodynam¬ 
ics can be treated as information processing under conser¬ 
vation laws, and so researchers began to experiment with 
other conserved quantities, like angular momentum [112— 
114], using the principle of maximum entropy to model 
thermal momentum baths. The state of those baths has 
again an exponential Gibbs form, with operators like L 
replacing H. The same type of monotones emerged, and 
similar behaviour was found for more general conserved 
quantities [ 11, 115]. 

9. Finite-size effects 

Another setting of practical interest is when we have 
access to a heat bath but may not draw arbitrary thermal 
subsystems from it. For instance, maybe we cannot cre¬ 
ate systems with a very large energy gap, or we can only 
thermalize a fixed number of qubits. In this case, the 
precision of state transformations is affected, as shown 
in [116], and we obtain effective measures of work cost 
that converge to the usual quantities in the limit of a 
large bath. The opposite limit, in which all resources 
are large heat baths, leads to the idea of heat engines 
(Example 2). 

10. Single-shot regime 

Some of the studies mentioned so far characterize the 
limit of many independent repetitions of physical exper¬ 


iments, and quantify things like the average work cost of 
transformations or conversion rates [86, 95]. The mono- 
tones found (like the von Neumann entropy and the usual 
free energy) are familiar from traditional thermodynam¬ 
ics, because this regime approximates the behaviour of 
large uncorrelated systems. As we move towards a ther¬ 
modynamic theory of individual quantum systems, it be¬ 
comes increasingly relevant to work in the single-shot 
regime. Some studies consider exact state transforma¬ 
tions [77, 78, 114], while others allow for a small error 
tolerance [79, 81, 82, 87, 88, 111, 115, 117, 118]. The 
monotones recovered correspond to operational entropy 
measures, like the smooth max-entropy (see Example 3) , 
and variations of a single-shot free energy that depend on 
the conservation laws of the setting; in general, they can 
be derived from quantum Renyi relative entropies [92] be¬ 
tween the initial state and an equilibrium state [83, 119]. 
Single-shot results converge asymptotically to the tradi¬ 
tional ones in the limit of many independent copies. The 
relation between single-shot and average regimes is stud¬ 
ied via fluctuation theorems in [120]. 


11. Definitions of work 

In classical thermodynamics, we can define work as 
some form of potential energy of an external device, 
which can be stored for later use. For instance, if a 
thermodynamic process results in the expansion of a gas 
against a piston, we can attach that piston to a weight, 
that is lifted as the gas expands. We count the gain in 
gravitational potential energy as work it is well-ordered 
energy that can later be converted into other forms, ac¬ 
cording to the needs of an agent. A critical aspect is that 
at this scale fluctuations are negligible, compared to the 
average energy gain. In the regime of small quantum sys¬ 
tems, this no longer holds, and it is not straightforward to 
find a good definition of work. Without a framework for 
resource theories of thermodynamics, a system for work 
storage is often left implicit. One option is to assume 
that we can perform any joint unitary operation U$b in 
a system S and a thermal bath B, and work is defined 
as the change in energy in the two systems manipulated, 
W := Tr (H sb Psb ) - Tr (H SB U S b Psb U\ b ), where 
Hsb is the (fixed) Hamiltonian of system and bath, and 
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Psb the initial state [116]. Another example, inheriting 
more directly from classical thermodynamics, assumes 
that we can change the Hamiltonian of S and bring it 
in contact with an implicit heat bath [121]; work at a 
time t is then defined as 

WW :=p t 'T r ( W (0^)- 

To study fluctuations around this average value, we con¬ 
sider work to be a random variable in the single-shot 
setting—this is explored in Section V. Note that in these 
examples work is not operationally motivated; rather it 
is defined as the change of energy that heat cannot ac¬ 
count for. Resource theories of thermodynamics, with 
their conservation laws, force us to consider an explicit 
system W for work storage. We act globally on S <S> W, 
and we can define work in terms of properties of the re¬ 
duced state of W. One proposal for the quantum equiva¬ 
lent of a weight that can be lifted, for the resource theory 
of thermal operations, is a harmonic oscillator, with a 
regular Hamiltonian H\y = J2n n e \ n ) ( n \- The energy 
gaps need to be sufficiently small to be compatible with 
the Hamiltonian of S; in the limit e —> 0 the Hamiltonian 
becomes Hw = J dx x |x) {x\ [87, 106]. Average work is 
defined as Tr (Hw Pvv al ) — Tr {Hw /o[v" ltlal ), and fluctua¬ 
tions can be studied directly in the final state of the work 
storage system, pw- This approach also allows us to ob¬ 
serve other effects, such as the build up of coherences 
in W, and of correlations between W and S. Another 
advantage is that we can adapt the storage system to 
other resource theories: for instance, we can have an an¬ 
gular momentum reservoir composed of many spins, and 
count work in terms of polarization of the reservoir [113]. 
These approaches are critically analysed in Ref. [122]; in 
particular, it is highlighted that they do not distinguish 
work from heat. For instance, thermalizing the work stor¬ 
age system may result in an increase of average energy, 
which is indiscriminately labelled as “average work”. In 
the same paper, an axiomatic approach to define work 
is proposed, based on concepts from resource theories 
and interactive proofs. There, work is seen as a figure of 
merit: a real function assigned to state transformations, 
W(p —>• a). Starting from a couple of assumptions, the 
authors derive properties of acceptable work functions 
W: for instance, that they can be written as the dif¬ 
ference between a monotone for initial and final state, 
W(p —» a) = g(p) — g(cr). The free energy is an example 
of such a valid work function. 


B. Generalizing resource theories 

Let us now abstract from particular resource theories, 
and think about their common features, and how we may 
generalize them. 


1. Starting from the pre-order 

As mentioned before, the set of allowed transforma¬ 
tions imposes a pre-order structure (5, <) on the state 
space S. One direction towards exploring the concept 
of resource theories could be to start precisely from 
such a pre-order structure. That was the approach of 
Caratheodory, then Giles and later Lieb and Yngvason, 
who pioneered the idea of resource theories for thermody¬ 
namics [77, 78, 126, 127]. In their work, the set of allowed 
transformations is implicitly assumed, but we work di¬ 
rectly with an abstract state space equipped with a pre¬ 
order relation. They were largely inspired by classical, 
macroscopic thermodynamics, as one may infer from the 
conditions imposed on the state space, but their results 
can be applied to thermodynamics of small quantum sys¬ 
tems [114]. Assuming that there exist minimal resources 
that can be scaled arbitrarily and act as ‘currency’, the 
authors obtain monotones for exact, single-shot state 
transformations. When applied to the pre-order relation 
on classical states that emerges from thermal operations, 
these monotones become single-shot versions of the free 
energy [114], 

2. Starting from the set of free resources 

In Ref. [119] general quantum resource theories are 
characterized based on the set of free resources of each 
theory. Assuming that the set of free states is well- 
behaved (for instance, that it is convex, and that the 
composition of two free states is still a free state), they 
show that the relative entropy between a resource and 
the set of free states is a monotone. This is because 
the relative entropy is contractive (non-increasing under 
quantum operations); the same result applies to any con¬ 
tractive metric. Finally, they find an expression for the 
asymptotic value of a resource in terms of this monotone: 
the conversion rate between two resources is given by the 
ratio between their asymptotic value. 

3. In category theory 

Ref. [1 ], and more recently Ref. [37] have generalized 
the framework of resource theories to objects known as 
symmetric monoidal categories. These can represent es¬ 
sentially any resource that can be composed (in the sense 
of combining copies of different resources, like tensoring 
states in quantum theory). The authors consider both 
physical states and processes as possible resources. After 
obtaining the pre-order structure from a set of allowed 
operations, resource theories can be classified according 
to several parameters. For instance, the authors iden¬ 
tify quantitative theories (where having more of a re¬ 
source helps, like for thermal operations) and qualitative 
ones (where it helps to have many different resources). 
They find expressions for asymptotic conversion rates in 
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Example 3: Landauer’s principle 


How much energy is needed to perform logical operations? What are the ultimate limits for heat dissipation of 
computers? These questions lie at the interface between thermodynamics and information theory, are of both 
foundational and practical interest. As Bennett realized, all computations can be decomposed into reversible 
operations followed by the erasure of a subsystem [123]. If we assume that the physical support of our computer 
is degenerate in energy, we recover the setting of noisy operations, in which unitaries are applied for free. That 
way, the thermodynamic cost of computation is simply the cost of erasure, which is defined as taking a system 
from its initial state p to a standard, predefined pure state |0) (like when we format a hard drive). Rolf 
Landauer first proposed that the work cost of erasing a completely unknown bit of information (think of a 
fully mixed qubit) in an environment of temperature T is fcgT In 2 [12]. That very same limit was also found 
for quantum systems, in the setting of thermal operations [116, 124], for the ideal case of an infinitely large 
heat bath and many operations; finite-size effects are analysed in Ref. [116]. 

Using Landauer’s principle as a building block, we can approach the more general question of erasing of a 
system that is not in a completely unknown state, but rather about which we have partial information. For 
example, imagine that we want to perform an algorithm in our quantum computer, and then erase a subsystem 
S (which could be a register or ancilla). The rest of our computer may be correlated with S, and therefore 
we can use it as a memory M, and use those correlations to optimize the erasure of S. In short, we want 
to take the initial state psm to |0) (0| s (g> pm, erasing S but not disturbing M. It was shown [ ’9, 82] that 
the optimal work cost of that transformation is approximately U^ lax (S|M) p fcBTln2, where e parametrizes our 
error tolerance and Hf nax (S\M) p is the smooth max entropy, a conditional entropy measure that measures our 
uncertainty about the state of S, given access to the memory M. It converges to the von Neumann entropy 
in the limit of many independent copies. In the special case where S and M are entangled, it may become 
negative—meaning that we may gain work in erasure, at the cost of correlations. Not incidentally, these results 
use quantum information processing techniques to compress the correlations between S and M before erasure; 
after all, ‘information is physical’ [125]. 


different regimes and, crucially, give varied examples of 
resource theories, within and beyond quantum theory, 
showing just how general this concept is. 

4. Resource theories of knowledge 

In Ref. [128], emphasis is given to the subjective knowl¬ 
edge of an observer. The framework introduced there al¬ 
lows us to embed macroscopic descriptions of reality into 
microscopic ones, which in turn lets us switch between 
different agents’ perspectives, and see how traditional 
large-scale thermodynamics can emerge from quantum 
resource theories like thermal operations. It also allows 
us to combine and relate different resource theories (like 
thermodynamics and LOCC), and to infer the structure 
of the state space (like the existence of subsystems or cor¬ 
relations) from modularity and commutativity of trans¬ 
formations. 


C. Outlook 

In the previous sections we identified several open 
problems. These can be grouped into two main direc¬ 
tions: 

• Quantumness: coherence, catalysis and clocks. It 
remains to find optimal coherent catalysts and 
clocks under realistic constraints (a generalization 


of Ref. [97]). This would give us a better under¬ 
standing of the thermodynamic power and limita¬ 
tions of coherent quantum states. It would also 
allow us to account for all costs involved in state 
transformations. 


• Identifying realistic conditions. We have been very 
good at defining sets of allowed transformations 
that are analogous to those of traditional thermo¬ 
dynamics, and recover the same monotones (like 
the free energy) in the limit of large, uncorrelated 
systems. The original spirit of thermodynamics, 
however, was to find transformations that were easy 
and cheap to implement for experimenters—for in¬ 
stance, those whose cost did not scale with the 
relevant parameters. In order to find meaningful 
resource theories for individual quantum systems, 
it is again imperative to turn to concrete experi¬ 
mental settings and try to identify easy and cheap 
transformations and resources. At this stage, it 
is not yet clear whether these will correspond to 
thermal operations, time-independent Hamiltoni¬ 
ans, or another model of quantum thermodynam¬ 
ics—in fact it is possible that they vary depend¬ 
ing on the experimental realization, from supercon¬ 
ducting qubits to ion traps. 





17 


IV. ENTANGLEMENT THEORY IN 
THERMODYNAMIC SETTINGS 

In the previous sections we have established how quan¬ 
tum information can be used to understand the very 
foundation of thermodynamics, from the emergence of 
thermal states to the resource theory of manipulating 
these with energy conserving unitaries. We have seen 
that phrasing thermodynamics as a resource theory can 
elucidate the meaning of thermodynamic quantities at 
the quantum scale, and how techniques originally devel¬ 
oped for a resource theory of communication can facil¬ 
itate this endeavour. The motivation behind this ap¬ 
proach is a very practical one: finding the ultimate lim¬ 
itations of achievable transformations under restrictions 
that follow from the nature of the investigated system 
that naturally limits the set of operations we can perform. 
As quantum information processing is becoming increas¬ 
ingly applied, we also need to think about fundamen¬ 
tal restrictions to quantum information itself, emerging 
from unavoidable thermodynamic considerations. There 
has thus been an increased interest in investigating sce¬ 
narios of quantum information processing where thermo¬ 
dynamic considerations cannot be ignored. From funda¬ 
mental limitations to the creation of QIP resources to 
their inherent work cost. In this section we try to give a 
brief overview over recent developments in this intersec¬ 
tion with a focus on the paradigmatic resource of QIP: 
entanglement. 


A. Correlations and entanglement under entropic 
restrictions 

Entanglement theory is in itself one of the most promi¬ 
nent examples of resource theories. Entanglement, a re¬ 
source behind almost all tasks in quantum information 
processing, is hard to create and once distributed can 
only decrease. Thus in entanglement theory classically 
correlated states come for free and local operations are 
considered cheap, which singles out entanglement as the 
resource to overcome such limitations. These limitations 
and resources are of course very different to the resources 
and tasks explored in the previous sections. A compre¬ 
hensive comparison between the principles behind these 
and more general resource theories is made in Ref. [129] 
and as examples of a more abstract treatment in Ref. [16]. 

Such resource theories are always designed to reflect 
specific physical settings, such as local operations and 
classical communication (LOCC) [130] as a natural con¬ 
straint for communication. It is therefore unavoid¬ 
able that when describing various physical circumstances 
these resource theories can be combined yielding hybrid 
theories. One natural example is the desire to process 
quantum information in a thermodynamic background. 
Ignoring limitations coming from available energies in a 
first step this leads to the task of producing resources for 
computation (such as entanglement or correlation) at a 


given entropy. Some of the first considerations in this di¬ 
rection were motivated by the prospect of using nuclear 
magnetic resonance (NMR) for quantum computation. 
Due to non-zero temperature, i.e. non-trivial restrictions 
on the entropy of the state, such systems would always 
be fairly close to the maximally mixed state. 

In this context the most natural question to ask, is 
whether a unitary transformation is capable of entan¬ 
gling a given input state. As a precursor to studying the 
possibility of entangling multipartite states, the complete 
solution for two qubits was found in Ref. [131] and later 
decent bounds on bipartite systems of arbitrary dimen¬ 
sion were presented in Ref. [132]. 

Another pathway was pursued by Refs. [133-135], 
where with NMR quantum computation in mind, vol¬ 
umes of separable states around the maximally mixed 
state were identified. These volumes imply that if any 
initial state is in close proximity of the maximally mixed 
state, there can be no chance of ever creating entangle¬ 
ment in such states, as the distance from the maximally 
mixed state is invariant under unitary transformations. 
Further improvements in terms of limiting temperatures 
were obtained in Ref. [136]. 

The question of whether a given state can be entan¬ 
gled under certain entropy restrictions clearly relies only 
on the eigenvalue spectrum of the considered state, as 
the best conceivable operation creating entanglement is a 
unitary one (which leaves eigenvalues unchanged). These 
questions were further pursued under the name of “sep¬ 
arability from spectrum in Refs. [137-139]. One of the 
main results important in the context of quantum ther¬ 
modynamics is the following: A state with eigenvalues 
Ai, ordered by size, i.e. {A^ > A^+i} can be entangled by 
an appropriate unitary if 

Ai > Ad_i + 2 \JXd-2^d • (10) 

More importantly, for 2 x m dimensional states, this con¬ 
dition is not only sufficient, but also necessary [139]. 

Moving beyond the mere presence of entanglement in 
the unitary orbit of input states, one encounters an in¬ 
trinsic difficulty of properly quantifying the entanglement 
created. There is a whole “zoo of entanglement mea¬ 
sures [143, 144] and only in the bipartite case there is a 
unique ’’currency” known, i.e. a paradigmatic resource 
state from which all other states can be created via LOCC 
(although recent progress has been made in the four qubit 
case, where it has been shown that after exclusion of a 
measure zero set, such a set of resource states can indeed 
be identified [145]). 

In any case one can at least study general correla¬ 
tions with a clear operational interpretation, such as 
the mutual information, which has been performed in 
Refs. [140, 146, 147]. In these papers the authors have, 
among other things, identified minimally and maximally 
correlated states in the unitary orbit of bipartite systems. 
It turns out that at least here the entropy poses only a 
rather trivial restriction and for any d-dimensional state 
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FIG. 4: Creating correlations between local thermal states 



Two systems, initially in local, product thermal states of temperature T, are correlated using a) any unitary 
[ 40], b) a unitary changing the average energy by at most A E [141] or c) a unitary on the systems and another 
bath at the same temperature T [142]. The local states of the two systems may heat up in the process: AT 
is infinite in a), a function of A E in b) and forced to be zero in c). In general, for setting b), the optimal 
conversion rate of average energy into mutual information is still unknown. 


p a mutual information of I p (A : B) = 21og 2 (d) — S(p) 
can be achieved via global unitary rotations. 

Exploiting these results Ref. [14: ] continued to study 
the generation of correlations and entanglement under 
entropic restrictions for multipartite systems. Inspired 
by the idea of thermal states as a free resource, the au¬ 
thors consider a multipartite system initially in a thermal 
state. They ask what is the highest temperature T ent at 
which entanglement can still be created, it scales with the 
dimension of the partitions and quantify the inherent cost 
in terms average energy change (see example (4) for an 
exemplary two qubit energy cost). By introducing con¬ 
crete protocols, i.e. unitary operations, the authors show 
that bipartite entanglement generation across all parti¬ 
tions of n-qubits is possible iff k B T/E < n/(21n(l-|-\/2)) 
and genuine multipartite entanglement across all parties 
can be created if k B T/E < n/(21n(n)) + 0(n/ ln(n) 2 ). 


B. Correlations and entanglement in a 
thermodynamic background 

In the context of thermodynamics the previous sub¬ 
section can be viewed as a very special case of operating 
on closed systems with an unlimited external energy sup¬ 
ply or a fully degenerate Hamiltonian. As elaborated in 
Section III of the review this does not encompass the 
whole potential of thermodynamic transformations. If 
the necessary correlating unitary does not conserve the 
total energy, we should account for the difference in av¬ 
erage energy between initial and final states. Taking into 
account also the average energy cost reveals an intrinsic 
work value of correlations and entanglement in general. 
This fundamental fact was first quantified in Ref. [ If]). 

Accounting for the average energy change in the uni¬ 
tary orbit of initial quantum states however still does 
not encompass the whole potential of thermodynamic re¬ 
source theories. Thermal operations on target states can 
also make use of a thermal bath at temperature T and 


thus can also reduce the entropy of the target system. 
Disregarding energy costs in this context of course yields 
the rather trivial result that any quantum information 
processing resource can be produced, simply by cooling 
the system (close) to the ground state and then perform¬ 
ing the adequate global unitary operation on it. Taking 
into account the free energy costs of correlating trans¬ 
formations, Ref. [142] has shown that every bit of corre¬ 
lation embodies an intrinsic work value proportional to 
the temperature of the system. For mutual information 
this yields a relation akin to Landauer’s principle for the 
work cost of creating correlations W cor , 

W cor > k B TI p (A : B), (11) 

and it implies a general free energy cost of entanglement 
that is bounded from above and below for the bipartite 
case in Ref. [142]. All previous considerations are illus¬ 
trated in Fig. 4. 

That extractable work can be stored in correlations is 
by no means a purely quantum phenomenon. Even clas¬ 
sical correlations can store work in situations where local 
work extraction is impossible. In Ref. [148] the quantum 
vs. classical capacity for storing extractable work purely 
in correlations was compared. For two qubits twice as 
much work can be stored in entangled correlations as the 
best possible separable (or even classical, which turns 
out to be the same) correlations admit (a fact that is 
also mentioned in Ref. [76] in a different setting). How¬ 
ever the difference between separably encoded work from 
correlations W sep and the maximally possible work in 
correlations W max scales as 

w^ = i - 0 ( 0 , ( 12 ) 

W max 

i.e. the quantum advantage vanishes in the thermody¬ 
namic limit of large systems. 

Concerning the extractable work from correlations one 
can also find seemingly contrary results if the figure of 
merit changes. The above considerations apply only if 
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Example 4: Entangling two qubits 


Creating entanglement from thermal states will always cost some energy. For the simplest case of entangling 
two qubits with energy gap E at zero temperature one can find a closed expression, e.g. for the concurrence, 
in terms of the invested average energy A E = W: 


C(W) = 



the target is an extraction of average energy or standard 
free energy, partially neglecting the details of the work 
distribution arising in the receiving system (detailed con¬ 
siderations of such work distribution fluctuations will be 
discussed in section (V)). One can just as well be inter¬ 
ested in a guaranteed amount of work. If that is the 
case one can arrive at more restrictions concerning work 
extraction as also recently demonstrated in [118]. Cu¬ 
riously in Ref. [149] it was shown, however, that these 
restrictions can be overcome by considering k initially 
uncorrelated catalysts that build up correlations in the 
process. In that context one can extract more determinis¬ 
tic work and can thus regard the stochastic independence 
of the input catalysts as a resource for work extraction, 
which is quite contrary to the case considered before and 
the thermodynamic limit. 

A different, but very related, setting exploring work 
gain from correlations is studied in the context of quan¬ 
tum feedback control. Here the task is rather to quantify 
the inevitable work cost arising from information gain in 
the process of a measurement. As in order to measure a 
system one needs to correlate with the system in question 
it follows intuitively that this scenario will also always in¬ 
duce work cost related to bipartite correlations between 
the system and the memory storing the information gain 
about the system. Here the work cost coming purely from 
correlations was quantified in Ref. [14 ], building upon 
older results on the inevitable work cost of quantum mea¬ 
surements [150-153] and Landauer’s principle. To model 
the necessary feedback control, the authors included a 
general model of a quantum memory upon which projec¬ 
tive measurements can be performed. The authors also 
studied the possible work gain from bipartite quantum 
states in this context. Denoting the state of the memory 
as pm the authors find an upper bound on the work gain 
(defined as the work extracted from both subsystems mi¬ 
nus the work cost of the measurements and subsequent 
erasure of the quantum memory) as 

W net < k B T(I p (A : B) - I p (A : B\p M )) - A F 0 (p). 

(13) 


C. Thermodynamics under locality restrictions 

In the previous subsections we have reviewed the 
prospect of creating quantum information processing re¬ 


sources in a thermodynamic background. The other ob¬ 
vious connection between the resource theories of entan¬ 
glement and thermodynamics is taking the converse ap¬ 
proach. Here one is interested in thermodynamic opera¬ 
tions under additional locality restrictions. 

In Ref. [15 ] the difference between the extractable 
work from bipartite quantum states in thermodynamics 
both with and without locality restrictions was studied. 
The resulting difference, called the work deficit , can be 
bounded via 

A = ma x[S(p A ),S(p B )] - S(pab) , (14) 

which for pure states coincides with entanglement of for¬ 
mation (or any other sensible choice of entanglement 
measure that all reduce to the marginal entropy in case 
of pure states). In the above equation it is assumed that 
bits which are sent down the communication channel are 
treated as classical in the sense that they are only de- 
pliased once, and not again in a second basis. This inter¬ 
play led to subsequent investigations into the thermody¬ 
namic nature of entanglement in Ref. [154], where analo¬ 
gies between irreversible operations in thermodynamics 
and bound entanglement were drawn, and to concrete 
physical scenarios satisfying this bound in Ref. [12 ]. 

D. Entanglement resources in thermodynamic 
tasks 

Apart from resource theory inspired questions, one 
might study the role of informational quantities through 
their inevitable appearance in thermodynamic operations 
at the quantum level. For instance the role of entan¬ 
gling operations and entanglement generation in extract¬ 
ing work from multiple copies of passive states , i.e. states 
where no local work extraction is possible [ 08, 109], has 
attracted some attention recently. The implied fact that 
global unitary operations are required to extract work 
indicates some form of non-local resource being involved 
in the process. 

In general passive states are always diagonal in the en¬ 
ergy eigenbasis [108, 109], which implies that one starts 
and ends the protocols with diagonal states. In these 
scenarios the individual batteries from which work is to 
be extracted are considered non-interacting, directly im¬ 
plying the separability of initial and final states in these 
protocols. Nonetheless the fact that local unitaries can 
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never extract any work from copies of passive states di¬ 
rectly implies that entangling unitaries enable work ex¬ 
traction from such states [155]. In that sense entangling 
power of unitaries can be seen as a resource for work ex¬ 
traction purposes (which in conventional thermodynamic 
resource is of course considered a free operation). 

In Ref. [156] the role of quantum resources in this con¬ 
text was further explored. While it is true that the abil¬ 
ity to perform entangling unitaries is required for this 
particular work extraction problem, this does not imply 
that any entanglement is ever generated in the process. 
In fact the whole procedure can dynamically be imple¬ 
mented without ever generating the slightest bit of en¬ 
tanglement [156], however the most direct transformation 
can considerably entangle the systems in the process. In 
Ref. [157] it was demonstrated that if the work per unit 
time (power) is considered with cyclic operations in mind 
then a quantum advantage for charging power can be 
achieved. 


E. Using thermodynamics to reveal quantumness 

That entanglement plays a special role in quantum 
many-body physics is a well established fact that has re¬ 
ceived adequate attention in numerous publications (see 
e.g. Ref. [158] and the extensive list of citations therein). 
In this topical review we want to at least mention a re¬ 
lated question that connects quantum thermodynamics 
directly with entanglement theory: The possibility to 
use thermodynamic observables to reveal an underlying 
entanglement present in the system. At zero tempera¬ 
ture it is already known that many natural interaction 
Hamiltonians have entangled ground states (in fact often 
many low energy eigenstates even of local Hamiltonians 
feature entanglement). This fact can be exploited to di¬ 
rectly use the energy of a system as an entanglement 
witness, even at non-zero temperatures [159]. Intuitively 
this can be understood through the fact that a low av¬ 
erage energy directly implies that the density matrix is 
close to the entangled ground state. If this distance is 
sufficiently small that can directly imply entanglement 
of the density matrix itself. The known results and open 
questions of this interplay including Refs. [160-167] are 
also discussed in the review Ref. [36]. Furthermore, other 
macroscopic thermodynamic quantities can also serve as 
entanglement witnesses through a similar intuition, such 
as e.g. the magnetic susceptibility [168] or the entropy 
[169], 

F. Outlook 

Resource theories always have their foundation in what 
we believe to be hard/impossible to implement and what 
resources allow us to overcome such limitations. As 
such they always only capture one specific aspect of the 
physical systems under investigation. The results out¬ 


lined in this section emphasise the fact that thermody¬ 
namic constraints have drastic consequences for process¬ 
ing quantum information and that locality constraints 
will change thermodynamic considerations at the quan¬ 
tum scale. One path to explore could now be a consis¬ 
tent resource theory that adaptively quantifies possible 
resources from different restrictions. This would further¬ 
more elucidate the exact role played by genuine quantum 
effects, such as entanglement, in thermodynamics. 


V. QUANTUM FLUCTUATION RELATIONS 
AND QUANTUM INFORMATION 

A. Introduction 

The phenomenological theory of thermodynamics suc¬ 
cessfully describes the equilibrium properties of macro¬ 
scopic systems ranging from refrigerators to black holes 
that is the domain of the large and many. By extrapo¬ 
lating backwards, from the domain of the ‘many’ to the 
‘few’, we venture further from equilibrium into a regime 
where both thermal and quantum fluctuations begin to 
dominate and correlations proliferate. One may then ask 
the question - what is an appropriate way to describe this 
blurry world which is dominated by deviations from the 
average behaviour? 

One way to describe the thermodynamics of small sys¬ 
tems where fluctuations cannot be ignored is by using 
the framework of stochastic thermodynamics [170]. In 
this approach the basic objects of traditional statistical 
mechanics such as work and heat are treated as stochastic 
random variables and hence characterised by probability 
distributions. Over the last 20 years various approaches 
have lead to sets of theorems and laws, beyond the lin¬ 
ear response regime, which have revitalised the already 
mature study of non-equilibrium statistical mechanics. 
Central to these efforts are the fluctuation relations that 
connect the non-equilibrium response of a system to its 
equilibrium properties. A wealth of results have been 
uncovered in both the classical and the quantum regimes 
and the interested reader is directed to the excellent re¬ 
views on the topics [21-23]. Here we focus on aspects of 
this approach that have been specifically influenced by 
concepts in quantum information, or show promise for 
symbiosis. We hope that by reviewing the existing con¬ 
tributions as well as suggesting possible research avenues, 
further cross fertilisation of the fields will occur. 

To begin with, it is useful to illustrate how the proba¬ 
bility distributions of a thermodynamic variable like work 
is defined. Consider a quantum system with a time- 
dependent Hamiltonian 77(A(i)), parametrized by the ex¬ 
ternally controlled work parameter A (t). The system is 
prepared in a thermal state by allowing it to equilibrate 
with a heat bath at inverse temperature /3 for a fixed 
value of the work parameter A (t < ti ) = Ai. The initial 


21 


state of the system is therefore the Gibbs state, 




e 

M Ai) 


At t = ti the system-reservoir coupling is removed and 
a fixed, reversible protocol is performed on the system 
taking the work parameter from its initial value A; to 
the final value Af at a later time t = if. The initial and 
final Hamiltonians are defined by their spectral decom¬ 
positions 


H{ Ai) =E„£n(Ai) \tpn) (0n| 

and H(\{) 10m) (0m | 

respectively, where | ij) n ) (|0 m )) is the nth (mth) eigen¬ 
state of the initial (final) Hamiltonian with eigenvalue 
E n (Xi),E m (X{). The protocol connecting the initial and 
final Hamiltonians generates the unitary evolution oper¬ 
ator U(tf : tf), which in general has the form 


= T_y exp 


-ij dt'H(X(t')) 


(15) 


where denotes the time-ordering operation. We 
stress here that, in this framework, one typically assumes 
that the system is initially prepared in a thermal state 
but after the unitary protocol the system is generally in 
a non-equilibrium state. 

The work performed (or extracted) on (or from) the 
system as a consequence of the protocol is defined by the 
outcomes of two projective energy measurements [171]. 
The first, at t = f,, projects onto the eigenbasis of the ini¬ 
tial Hamiltonian H( Af), with the system in thermal equi¬ 
librium. The system then evolves under the unitary op¬ 
erator U(tf 7 ti) before a second projective measurement is 
made onto the eigenbasis of the final Hamiltonian H (Af) 
at t = tf. The joint probability of obtaining the outcome 
E n { Af) for the initial measurement followed by E m ( Af) 
for the final one is easily shown to be 


-/3B„(Ai) 

p(n,m)= 2 |(0m|f7(Mi)|0n)| • (16) 


Accordingly, the quantum work distribution is defined as 


P F (W) = y^p(n,m) S(W - [E m ( A f ) - E n (Xi )]). 

n,m 

(17) 


where <5 is the Dirac delta function. For reasons which 
will become clear shortly we use the subscript F to de¬ 
note ‘forward’ protocol. Physically, Eq. 17 states that 
the work distribution consists of the discrete number of 
allowed values for the work (E m ( Af) — -E„(Ai)) weighted 
by the probability p(n, m) of measuring that value in a 
given realisation of the experiment. The quantum work 
distribution therefore encodes fluctuations in the mea¬ 
sured work arising from thermal statistics (first measure¬ 
ment) and from quantum measurement statistics (second 
measurement). 


In order to understand what is meant by a fluctua¬ 
tion theorem, we introduce a backward process which 
is the time reversed protocol of the forward one previ¬ 
ously defined. Now P B (W) is the work distribution cor¬ 
responding to the backward process , in which the system 
is prepared in the Gibbs state of the final Hamiltonian 
H( Af) at t = 0 and subjected to the time-reversed proto¬ 
col that generates the evolution QU(t{, ti)© 1 ', where 0 is 
the anti-unitary time-reversal operator. It turns out that 
the following theorem holds, the Tasaki-Crooks relation 
[172, 173], 


Pf{W) _ p(W-AF) 

Pb(-W) 


(18) 


which shows that, for any closed quantum system under¬ 
going an arbitrary non-equilibrium transformation, the 
fluctuations in work are related to the equilibrium free 
energy difference for the corresponding isothermal pro¬ 
cess between the equilibrium states r(Af) and r(A/), 


AF P 


In 


f Zp(Xj) \ 

U(A f) ;- 


(19) 


This relationship is further emphasized by a corollary 
to Eq. 18 known as the Jarzynski equality [174], 

J dWP F (W)e~ pw = (e~P w ) = e~^ AF (20) 


which states that AF (of the corresponding isother¬ 
mal process) can be extracted from by measuring the 
exponentiated work. A straightforward application of 
Jensen’s inequality for convex functions allows one the 
retrieve the expected expression ( W ) > AF. The aver¬ 
age energetic deviation of a non-equilibrium process from 
the equivalent reversible isothermal process is known as 
dissipated work 


(W) diss = (W) - AF. (21) 


Due to the Jarzynski equality this quantity is positive, 
( W)diss > 0. This can be also directly seen from the 
Crooks relation, taking the logarithm of both sides of 
the equality in Eq. 18 and integrating over the forward 
distribution we find 


(E) = (3{(W) - AF) = K(P F (W)\\P B (-W)) (22) 

where K is the classical Kullback Leibler divergence 
and we have introduced the average irreversible entropy 
change (E) corresponding to the dissipated work. Phys¬ 
ically the irreversible entropy change, in this context, 
would be the internal entropy generated due to the non¬ 
equilibrium process which would manifest itself as an ad¬ 
ditional source of heat if an ideal thermal bath would 
be reconnected to the system at the end of the protocol. 
In Ref. [175] it was shown that the irreversible entropy 
change can also be expressed in terms of a quantum rel¬ 
ative entropy 

(E) = D(ct||t(A / ,/3)) (23) 

where a = U(tf, t;)r(Aj, fi)U\tf, ti) is the out of equi¬ 
librium state at the end of the protocol. This is fully 
consistent with the open system treatment in [176]. 
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B. Phase estimation schemes for extraction of 
quantum work and heat statistics 

Surprisingly, perhaps one of the most important contri¬ 
butions that ideas from quantum information have made 
to this field in statistical mechanics is the experimental 
acquisition of statistics of work. In the classical setting 
considerable progress has been made in the experimen¬ 
tal extraction of the relevant stochastic thermodynamic 
distributions to explore and verify the fluctuation theo¬ 
rems [23]. Up until very recently, no such experimental 
progress had been made for quantum systems. A central 
issue is the problem of building the quantum work distri¬ 
bution as it requires to make reliable projective energy 
measurements on to the instantaneous energy eigenbasis 
of an evolving quantum system [22, 171]. It was proposed 
in Ref. [177] that these measurements could be reliably 
performed on a single trapped ion, an experiment that 
was recently performed [178]. 

Alternatives to the projective method have been pro¬ 
posed [179, 180], based on phase estimation schemes , well 
known in quantum information and quantum optics [181]. 
In these schemes, we couple our system to an ancillary 
system, and perform tomography on that system. The 
spirit is very similar to the DQC1 algorithm put forward 
in Ref. [182]. The characteristic function of the work 
probability distribution (Eq. 24) can be obtained from 
the ancilla, and the work statistics are then extracted by 
Fourier transform. The characteristic function is defined 
as 

XF (u) = j dW e iuW P F (W). (24) 

The proposals to measure the characteristic function were 
first tested in the laboratory only quite recently in a Liq¬ 
uid state NMR setup [183]. This experiment is the first 
demonstration of the work fluctuation theorems and ex¬ 
traction of work quantum statistics, and is expected to 
inspire a new generation of experiments at the quantum 
level. Another interesting extension of these schemes is to 
go beyond the closed system paradigm and to study open 
system dynamics at and beyond the weak coupling limit. 
The first extensions have been proposed in Refs. [184- 
186]. In Ref. [187] the proposal outlined in Ref. [186] 
to measure the statistics of dissipated heat was imple¬ 
mented in order to perform a study of the information 
to energy conversion in basic quantum logic gates at the 
fundamental Landauer Limit. 

Another interesting suggestion made to access the 
quantum work statistics is the use the concept of a ‘pos¬ 
itive operator valued measure’, or POVM [188], a well- 
known concept within quantum information and quan¬ 
tum optics. A POVM is the most general way to describe 
a measurement in a quantum system, with the advantage 
that it can always be seen as a projective measurement on 
an enlarged system. In this work the authors show that 
by introducing an appropriate ancilla that the POVM 
description allows the work distribution to be efficiently 


sampled with just a single measurement in time. In this 
work it was suggested that the algorithm proposed could 
be used, in combination with the fluctuation theorems, 
to estimate the free energy of quantum states on a quan¬ 
tum computer. The scheme was recently extended and 
developed in Ref. [189] along with a promising implemen¬ 
tation using ultra-cold atoms. This would be a promising 
avenue to explore work statistics in a many-body physics 
setting where the statistics of work can be shown to have 
universal behaviour at critical points [190]. 

C. Fluctuation relations with feedback, 
measurement and CPTP maps 

The relationship between thermodynamics and the in¬ 
formation processing is almost as old as thermodynamics 
itself and is no where more dramatically manifested than 
by Maxwell’s demon [2-6] . One way of understanding the 
demon paradox is by viewing the demon as performing 
feedback control on the thermodynamic system. In this 
case the framework for stochastic thermodynamics and 
the fluctuation theorems needs to be expanded. Build¬ 
ing upon previous work [152, 153], Sagawa and Ueda 
have generalised the Jarzynski equality to incorporate 
the feedback mechanism [191, 192] for classical systems. 
This theoretical breakthrough allowed for an experimen¬ 
tal demonstration of information to energy conversion in 
a system by means of of non-equilibrium feedback of a 
Brownian particle [193]. These feedback based fluctua¬ 
tion theorems were further modified to incorporate both 
initial and final correlations [194], These works, in par¬ 
ticular, highlight the pivotal role played by mutual infor¬ 
mation in non-equilibrium thermodynamics [6]. 

The Sagawa-Ueda relations were generalized to quan¬ 
tum systems in Ref. [195]. For reasons of pedagogy we 
will follow this approach here. In the work of Morikuni 
and Tasaki an isolated quantum system is considered 
where an external agent has control of the Hamiltonian 
parameters. The system is initialised in a canonical state, 
r(/3), and an initial projective measurement of the energy 
is made whose outcome is E The Hamiltonian is then 
changed via a fixed protocol and evolves according to the 
unitary operator U. In the next stage a projective mea¬ 
surement is performed with outcomes j = 1,... ,n and 
described by a set of projection operators Hi,...,n„. 
Now the time evolution is conditioned on the outcome 
j so the Hamiltonian is changed according to these out¬ 
comes. This is the feedback control stage. Finally, one 
makes a projective measurement of the energy of the fi¬ 
nal Hamiltonian with outcome E J k . In this setting it is 
shown that 

{e P(W- AF)> = 7) (25) 

where W = Wij t k = Ef — E k is the work and A F is the 
free energy difference between the initial state and the 
canonical state corresponding to the final value of the 
Hamiltonian W. We see that in this feedback controlled 
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scenario a new term enters on the right hand side. A 
straightforward calculation shows that this term evalu¬ 
ates as 7 = JATr[II jU^r{P)UjTlj\. This 7 quantity is 
shown in Refs. [191, 195] to be related to the efficiency of 
the demon in making use of the information it acquires 
during the feedback process. When it becomes less than 
one it provides an example of a failed demon who did 
not make a good use of the information acquired. On 
the other hand it can become larger than one indicating 
that the feedback is working efficiently. Another rela¬ 
tion discovered by Sagawa and Ueda and quantized by 
Morkikuni and Taskaki concerns almost the same proto¬ 
col as just explained only now classical errors are made in 
the intermediate measurement stage. Again let the inter¬ 
mediate measurement be described by IIi,...,n n which 
yield the result j but the controller misinterprets the re¬ 
sult as j' with a certain probability. In this framework 
another generalised fluctuation theorem can be derived, 

{SW-AF-I )) = J (26) 

where / is the mutual information between the set of 
measurement outcomes the demon actually records and 
what is the true result of the projection. These feedback 
fluctuation theorems for quantum systems were further 
generalised to the situation when a memory system is 
explicitly accounted for in Ref. [14 ] and shed light on 
the amount of thermodynamic work which can be gained 
from entanglement. In addition to feedback, fluctuation 
theorems were investigated under continuous monitoring 
[196, 197] and analysed for general measurements [198, 
199]. 

A recent series of papers have analysed fluctuation¬ 
like relations from the operational viewpoint employing 
the full machinery of trace-preserving completely posi¬ 
tive maps. In Ref. [200] the formalism is used to give 
an alternative derivation of the Holevo bound [201]. In 
Ref. [202] an information-theoretical Jarzynski equality 
was derived. It was found that fluctuation relations can 
be derived if the map generated by the open dynamics 
obeys the unital condition. This has been connected to 
the breakdown of micro-reversibility for non-unital quan¬ 
tum channels [203-206]. In Ref. [207] the authors anal¬ 
ysed the statistics of heat dissipated in a general proto¬ 
col and found that the approach can be used to derive a 
lower bound on the heat dissipated for non-unital chan¬ 
nels. Recently this bound has been used to investigate 
the connection with the build up of multipartite correla¬ 
tions in collisional models [208]. 

D. Entropy production, relative entropy and 
correlations 

With the surge of interest in the thermodynamics of 
quantum systems and the development of quantum fluc¬ 
tuation relations, research has been directed to micro¬ 
scopic expressions for entropy production. In formulat¬ 
ing thermodynamics for non-equilibrium quantum sys¬ 


tems, the relative entropy plays a central role [192]. As 
first pointed out in Ref. [209] this is due to its close re¬ 
lationship with the free energy of a quantum state. The 
relative entropy also plays a central role in quantum in¬ 
formation theory, in particular, in the geometric picture 
of entanglement and general quantum and classical corre¬ 
lations [210, 211]. In the non-equilibrium formulation of 
thermodynamics [ 22 ] it is omnipresent for the description 
of irreversible entropy production in both closed [175] 
and open driven quantum systems [212] (see also [213]). 
One may then wonder if there exists a relationship be¬ 
tween the entropy produced by operations that generate 
or delete correlations in a quantum state and the mea¬ 
sures for correlations in that state? Given the youthful 
nature of the field the question is largely unanswered but 
some progress in this direction has been made. 

The relationship between the relative entropy of en¬ 
tanglement and the dissipated work was first proposed 
as an entanglement witness in Ref. [214]. Going beyond 
the geometric approach a functional relationship between 
the entanglement generated in a chain of oscillators and 
the work dissipated was explored in Ref. [215] and also 
later for more general quantum correlations [216]. In 
an open systems framework it was shown that the ir¬ 
reversible entropy production maybe attributed to the 
total correlations between the system and the reservoir 
[217] (we note that this derivation is entirely analogous 
to the formulation of the Landauer principle put forward 
by in Ref. [116]). The exchange fluctuation relation and 
the consequences for correlated quantum systems were 
studied in Ref. [218]. 


E. Outlook 

As fluctuation theorems are exact results, valid for ar¬ 
bitrary non-equilibrium dynamics, they are currently be¬ 
ing used to understand the non-linear transport of en¬ 
ergy, heat and even information in quantum technologies. 
This is a relatively new research avenue and the appli¬ 
cations of quantum fluctuation theorems in other fields 
such as condensed matter physics, quantum optics and 
quantum information theory are in their infancy. Ulti¬ 
mately, the hope would be that they provide a unifying 
framework to understand the relationship between infor¬ 
mation and energy in non-equilibrium quantum systems. 
Ultimately one would like to form a picture of informa¬ 
tion thermodynamics of quantum systems under general 
non-equilibrium conditions. 

As we have seen above, quantum phase estima¬ 
tion, a central protocol in quantum information theory, 
has been applied successfully to extract work statistics 
from a small non-equilibrium quantum system and per¬ 
haps other such unexpected interdisciplinary links will 
emerge. For example one wonders if existing experimen¬ 
tal schemes could be modified to deal with situations 
dealing with non-passive initial states so as to study max¬ 
imal work extraction problems and also to extend to more 
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complicated many-body and open system scenarios. 

In Refs. [120, 219, 220] the first steps towards unifica¬ 
tion of the work statistics and fluctuation theorems ap¬ 
proach to thermodynamics and the single shot statistical 
mechanics approaches mentioned have been taken (see 
Sec. III). We are confident that other links will emerge 
between various approaches in the not so distant future. 


VI. QUANTUM THERMAL MACHINES 

In this final brief section of the review we end by con¬ 
sidering the area of quantum thermodynamics concerning 
quantum thermal machines, that is quantum versions of 
heat engines or refrigerators. We shall overview the ex¬ 
tent to which quantum entanglement and correlations are 
relevant to their operation. 

Whereas in almost all of the above the situation com¬ 
prised of only one thermal bath and systems in contact 
with it, in this section our interest is in situations involv¬ 
ing two (or more) thermal baths. Now, there are two 
regimes which one can focus on: the primary one is usu¬ 
ally the cyclic behaviour of systems interacting with the 
baths, or alternatively the steady state behaviour that is 
characterised by the currents of heat or work that can be 
maintained in the long time limit. The second regime is 
the transient one, and how the system reaches stationar- 
ity. 

One way to think of the present situation is that the 
second thermal bath is the system out of equilibrium with 
respect to the first bath, and the goal is to produce re¬ 
sources (work, or a steady state current out of a cold 
bath) at optimal rates. From this perspective, the quan¬ 
tum machine plays the role of the ‘bridge’ or the ‘media¬ 
tor’ which facilitates the operation of the larger thermal 
machine. 

The history of quantum thermal machines is a long 
one, going back to the sixties with the invention of the 
maser, which can be seen as a heat engine [ 221 ], and 
received much attention over the following decades. A 
complete overview of the literature in this direction is far 
beyond the scope of the present review; however excel¬ 
lent recent overviews can be found in Refs. [31-33]. In the 
present context, one important message from this body 
of work is that thermal machines comprised of as little as 
a single qutrit (3 level system), or of 2 or 3 qubits, can 
be constructed, that moreover can approach Carnot effi¬ 
ciency (the maximal possible efficiency of any machine). 
It is thus plausible that they may ultimately become im¬ 
portant from the perspective of nanotechnology and im¬ 
plementations of quantum information processes devices. 
As such a full understanding of their quantum behaviour, 
including the correlations they can build up, is impor¬ 
tant. Here we review specifically those studies concerned 
with the role of entanglement and quantum coherence in 
the functioning of such small quantum thermal machines, 
both at the level of the machine, as well as in the bath, 
if pre-processing operations are allowed. We also look at 


the role of coherence in the transient behaviour when the 
refrigerator is first switched on. We review a recent pro¬ 
posal for a witness that quantum machines are provably 
outperforming their classical counterparts. Finally, we 
look at the idea of using thermal machines as a means of 
entanglement generation (switching the focus away from 
the traditions resources of work or heat currents). 

A related idea is that of algorithmic cooling , which 
we summarise in Example 5, and which was recently re¬ 
viewed in [ 222 ]. 


A. Absorption refrigerators 

The first machine we shall look at a quantum model 
of an absorption refrigerator, a refrigerator which is not 
run by a supply of external work (which is the situa¬ 
tion most customarily considered), but rather run by a 
source of heat. An absorption refrigerator is thus a device 
connected to three thermal reservoirs; a ‘cold’ reservoir 
at temperature /3c from which heat will be extracted; 
a ‘hot’ reservoir at inverse temperature /3 h, which pro¬ 
vides the supply of energy into the machine; and finally 
a ‘room temperature’ reservoir at temperature /3 r into 
which heat (and entropy) will be discarded. The goal is 
to cool down the cold reservoir (i.e. extract heat from 
it). 

There are a number of different figures of merit that 
one can consider to quantify the performance of the ma¬ 
chine. The most commonly considered is the coefficient 
of performance COP = Qc/Qn , where Q c and Qr are 
respectively the heat currents flowing out of the cold the 
hot reservoirs (the COP is the analogous quantity to the 
efficiency for an absorption refrigerator; since the COP 
can be larger than 1 it cannot be thought of directly as 
an efficiency). The famous result of Carnot [223] , a state¬ 
ment of the second law of thermodynamics, is that the 
efficiency (or COP) of all thermal machines is bounded 
as a function of the reservoir temperatures. In particular, 
for the specific case of an absorption refrigerator we have 
COP < (/3 r — /3h)/(/3c — /3r). Other relevant figures of 
merit are the power Qc (i.e. neglecting how efficient the 
process is), the COP when running at maximal power, 
and the minimal attainable stationary temperature /3q 
for a cold object in contact with the bath. 

Below we give a brief outline of the model under 
consideration, full details of which can be found in 
Refs. [224, 221 ]. Consider three qubits, each one in ther¬ 
mal contact with one of the three thermal baths, with 
local Hamiltonians Hi = Ei |1) (1|, for i = C, R, H chosen 
such that Er = Eq + Ar to ensure that the system has 
a degenerate subspace of energy Er formed by the states 
|010) and 1101) (where we use the order C-R-H for the 
three qubits). In this subspace the interaction Hamilto¬ 
nian H lnt = <?(|010) (1011 + 1101) (010|) is placed, which 
mediates the transfer of energy. A schematic representa¬ 
tion of this fridge can be found in Fig. 5. 
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Example 5: Algorithmic cooling 


Consider a collection of n qubits, all at inverse temperature /3, with corresponding populations in the ground 
and excited states p and (1 —p) respectively. The goal of algorithmic cooling is to bring m qubits to the ground 
state by an arbitrary unitary transformation. A fundamental upper bound can be placed on m, purely by 
entropic considerations. The initial entropy is nS(r(l 3) = nH(p ), where H(p) = — p\og 2 p — (1 — p) log 2 (l — p) 
is the binary Shannon entropy. Since unitary transformations do not change the entropy, this easily leads to 
the upper bound on m, 

m<n(l — H(p)) (27) 

which would be achieved if the remaining n — m qubits are all left at infinite temperature (maximally mixed 
state) with entropy 5(r(0)) = 1. In [226] it was shown that as n tends to infinity this fundamental limit can 
be approached using an algorithm which uses 0{n log 2 n) unitary gate operations. It was later realised that 
given access to an external bath this limit can be surpassed: the qubits which end this protocol at infinite 
temperature can be ’refreshed’ to temperature /3 and the protocol can be run again on the remaining (n — m) 
qubits, for example [227]. This is referred to as heat-bath algorithmic cooling. 

In order to understand the basic principle, one can focus instead on 3 qubits and assume that the first is the 
one which is to be cooled down (now not to zero temperature, but any colder temperature). Let us consider 
the populations of the two states 1100) and |011), which are p 2 (l— p) and p(l—p) 2 respectively. The state|100), 
in which qubit one is excited (and therefore ‘hot’) has more population than the state 1011), where qubit one 
is in the ground state (and therefore ‘cold’). Thus, by swapping the population of these two states the first 
qubit is cooled down. Indeed, after the application of such a unitary, the final population p' in the ground 
state of the first qubit is 

P 1 =P+ (2p- l)p(l -p) (28) 

which is greater than p whenever (2 p — 1) > 0, i.e. whenever the first qubit was at a positive temperature. 
Finally, a unitary which implements |011) O 1100) whilst leaving all other energy eigenstates the same can 
easily be constructed from the CNOT and Toffoli gates as 
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A recent review giving many more details about algorithmic cooling can be found in Ref. [222] 


1. Stationary behaviour 

Assuming the weak coupling regime between the qubits 
and the baths, the dynamics can be modelled using a 
time-independent Lindblad Master equation p = C(p ) 
(with C the Linbladian, i.e. the most general generator 
of time-homogeneous, Markovian dynamics). The sta¬ 
tionary solution p st , satisfying C(p st ) = 0, can be shown 
to correspond to an absorption refrigerator if the param¬ 
eters are chosen appropriately, i.e. such that /3 q > /3c, 
where /3 q is the stationary inverse temperature of the 
cold qubit. 

From the point of view of quantum information, the ba¬ 
sic questions about this steady state are (i) whether quan¬ 
tum correlations (for example entanglement) are present 
in the stationary state, and (ii) if yes, whether they are 
important for the operation, or merely a by-product of 
quantum evolution. These questions were addressed in 
Refs. [224, 225], 


In Ref. [225] quantum correlations in the form of dis¬ 
cord were studied. The quantum discord T>(AB) p := 
I(A : B) p — I(A : B) ai with a the state after a minimally 
disturbing measurement on Bob, is a form of quantum 
correlation weaker than entanglement [228, 229]. The 
authors studied quantum discord between numerous in¬ 
equivalent partitions of the system. The most interesting 
results were obtained when the discord is calculated be¬ 
tween the cold qubit (the qubit which is being cooled) 
and the relevant subspace of the two remaining qubits 
(that singled out by the interaction Hamiltonian Hi nt ). 
They found that discord is always present, but they found 
no relationship between the amount of discord present 
and the rate at which heat was extracted from the cold 
bath. Specifically, to obtain this result they studied the 
behaviour of discord as a function of the energy spacing 
Eq of the cold qubit. Whilst both quantities typically 
exhibited local maxima as Ec was varied, these maxima 
failed to coincide. 

In Ref. [224] the focus was instead on the entanglement 
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FIG. 5: Three-qubit fridge 



Schematic diagram of a three qubit autonomous refrigerator (inside circle), coupled to three thermal reservoirs. 
The interaction Hamiltonian is represented by the green and orange arrows. 


maintained in the steady state. First, if the machine is 
operating close to the maximal Carnot limit then the 
state is necessarily fully separable, i.e. a convex combi¬ 
nation of product states of the three qubits. Conversely, 
operating far from this regime every type of multipar¬ 
tite entanglement can be found in the stationary state. 
In particular, there are regimes where entanglement is 
generated across any fixed bipartition, and even genuine 
multipartite entanglement can be found, demonstrating 
that the state has no biseparable decomposition. Here it 
must be stressed that the amount of entanglement found 
was small, but that this should be expected due to the 
weak inter-qubit coupling. 


Finally, it was also shown that there appears to be 
a link between the amount of entanglement generated 
in the partition R|CH and the so-called cooling advan¬ 
tage that entangled machines have compared to separa¬ 
ble ones. In particular, the cooling advantage was defined 
as the difference between the minimal possible tempera¬ 
tures that could be achieved with either separable or en¬ 
tangled refrigerators. More precisely, by optimising the 
stationary temperature /3q of the cold qubit, varying the 
Hamiltonian of the machine qubits and their couplings to 
baths (at fixed temperatures). It was shown that arbi¬ 
trary machines (i.e. ones allowed to be entangled) could 
outperform ones which were additionally constrained to 
be separable. Moreover, the advantage was found to be 
a function only of the amount of entanglement generated 
across the R|CH partition. One point of interest is that 
this is the bipartition of energy entering vs. energy leav¬ 
ing the machine, thus suggesting a connection between 
the transport properties of the machine and the entan¬ 
glement. 


2. Transient behaviour 

Instead of looking at the steady state behaviour, one 
may also consider the transient behaviour. Such ques¬ 
tions are relevant when one is interested in running a 
small number of cooling cycles in order to cool down 
the system as fast as possible. Alternatively, if one is 
thinking of initialising a system for some other use, the 
transient regime might also be of interest for quicker ini¬ 
tialisation. Intuitively, since the evolution between the 
qubits is coherent, one might expect the local popula¬ 
tions to undergo Rabi oscillations, and hence by running 
for precise times lower temperatures may be achievable 
in a transient regime (as the qubits continuous cool down 
and heat up). 

This is precisely what was shown in Ref. [230, 231]. 
More precisely, in Ref. [231] the authors study the Marko¬ 
vian dynamics with weak inter-qubit coupling g (rela¬ 
tive to the relaxation rates, as in the above subsection), 
while in Ref. [230] the authors considered additionally 
Markovian dynamics with strong inter-qubit coupling, 
and band-limited non-Markovian baths (modelled with 
a one-qubit memory for each machine qubit). Taking as 
the natural initial state the product state with qubit to 
be initially at the same temperature of the bath, both 
numerically study the transient behaviour of the tem¬ 
perature of the cold qubit as the system approaches sta- 
tionarity. While in the weak interaction case no Rabi 
oscillations are observed (since the system is effectively 
over-damped), in the strong-interaction case Rabi oscilla¬ 
tions indeed take place, with period approximately 2n/g. 
This demonstrates that coherent oscillations offer an ad¬ 
vantage for cooling. A more complicated behaviour due 
to memory effects is also observed in the non-Markovian 
case in [230] , but nevertheless the system can be seen to 
pass through much colder temperatures during its tran¬ 
sient behaviour. In Ref. [231] it was also shown that 
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if the couplings are chosen appropriately, (in particular 
such that the weakest coupling is to the hot reservoir), 
then the system can quickly remains for a long time in 
a temperature below the stationary temperature, in par¬ 
ticular without oscillating above it. This demonstrates a 
particular stable regime for the preparation of the system 
at temperatures below its stationary temperature. 

In order to explore more the advantage offered by co¬ 
herence, Ref. [230] also considered varying the initial 
state, by altering the coherence in the subspace where the 
Hamiltonian operates. Interestingly, with only a small 
amount of initial coherence even when considering case 
(a) of weak-interaction dynamics, oscillations in the tem¬ 
perature are seen, again allowing for cooling below the 
stationary temperature. In the other two cases, the mag¬ 
nitude of the oscillations is also seen to increase (i.e. the 
system achieves lower temperatures transiently), demon¬ 
strating an advantage in all situations. 

Finally, in Ref. [231] the amount of entanglement that 
is generated in the transient regime was also studied. Fo¬ 
cusing on either genuine multipartite entanglement, or 
entanglement across the partition R|CH, i.e. the one 
corresponding to energy-in vs. energy out (as studied 
in Ref. [224]), considerably more entanglement can be 
generated in the transient regime. 


B. Reservoir engineering 

As we have seen in previous sections of the review, 
thermals states are naturally considered as a free resource 
which can be utilised and manipulated. Likewise, the 
ubiquity of thermal machines is that having access to 
two large thermal reservoirs can also be considered as 
something essentially free, and thermal machines con¬ 
sider ways of utilising these resources. 

One interesting avenue is to consider that any trans¬ 
formation of a thermal reservoir which can ‘easily’ be 
carried out can also be considered to be free, as an ideal¬ 
isation, and this motivates the idea of considering ther¬ 
mal machines which run between engineered reservoirs , 
assuming that the engineering was an easy to perform 
transformation. In the present context, when one has 
sufficient control over (part of) the reservoir, then the 
engineering can be at the quantum level. Here again we 
are interested specifically in the role that quantum cor¬ 
relations engineered in the bath have on the functioning 
of quantum thermal machines. 

In Refs. [232, 233] reservoir engineering in the form of 
squeezing is considered, since squeezing is relatively easy 
to carry out, and is furthermore known to offer quan¬ 
tum advantages in other contexts in quantum informa¬ 
tion. That is the reservoir, instead of consisting of a large 
collection of modes in thermal states at inverse tempera¬ 
ture /3 h , are in fact squeezed thermal states (at the same 
temperature). More precisely, the squeezing operator is 
I/ sq = exp((ra 2 — r*a^ )/2) with a and o' the annihilation 
and creation operators respectively, and the squeezed 


thermal state (of a given mode, i.e. a harmonic oscil¬ 
lator) is U sq T(B)Uj q . Whereas normally the variances of 
the quadratures (x = (a + a')/2 and p = (a — a')/2i) 
are symmetric, the squeezed modes become asymmet¬ 
ric, with one the former amplified by the factor e r , and 
the latter shrunk by e~ r . The important point is that a 
system placed in thermal contact with such a squeezed 
reservoir will not thermalize towards a thermal state at 
P, but rather to a squeezed thermal state, which has the 
same average number of photons as a thermal state at 
temperature j3{r) < p. That is, in terms of average num¬ 
ber of photons, a squeezed thermal state appears ‘hotter’ 
than a thermal bath. 

Starting first with Ref. [233], a model of an absorption 
refrigerators is considered, identical to the one outlined 
in the previous section. Here, in accordance with the 
above, in the weak coupling regime the effect of the reser¬ 
voir engineering amounts to modifying the Linbladian C , 
such that the term corresponding to the hot reservoir 
£h transforms to Cn(r), where this now generates dis¬ 
sipation towards the squeezed thermal state at Pr {r) ■ 
They show that maximal COP that the refrigerator can 
approach becomes 


V(r) 


Pr - Pn{r) 
Pc ~ Pr 


> Vc 


Pr - Pr 
P c - Pr 


(29) 


That is, the COP overcomes the Carnot limit that 
bounds the COP of any absorption refrigerator operating 
between baths at Pc, Pr and Pr, if reservoir engineer¬ 
ing is not carried out. Thus if reservoir engineering is 
more readily available than a hotter ‘hot’ bath, then this 
approach clearly provides an advantage in terms of COP. 

In Ref. [232] a different model was considered, this time 
a quantum heat engine operating a quantum Otto cycle, 
a time dependent cycle, comprising two expansion stages 
(changing the Hamiltonian of the system) and two ther- 
malization stages. This system considered comprised of a 
single harmonic oscillator, with initial spacing E \. While 
uncoupled to any environment, the first stage is an ex¬ 
pansion, whereby E\ —> E 2 > E\, i.e. the Hamiltonian is 
changed in time. In the second stage the system is then 
placed in contact with a squeezed hot reservoir (this is the 
stage which differs from a standard Otto cycle, where an 
unsqueezed hot reservoir is used). After disconnection, 
the third stage is a compression stage, bringing the spac¬ 
ing back to from E 2 to Ey. Finally, the system is placed 
in contact with a cold (unsqueezed) reservoir, in order to 
thermalize at the cold temperature. This cycle is sum¬ 
marised in Fig. 6. The authors perform an analysis of 
the system and similarly show that the maximum effi¬ 
ciency of the engine exceeds the Carnot efficiency (of the 
Otto cycle, ij = Pr/Pc)- Moreover, if one considers the 
efficiency at maximum power, then this can also be sur¬ 
passed, and as the squeezing parameter becomes large, 
the efficiency at maximum power approaches unity. 

Finally, we stress that these results do not constitute 
a violation of the second law, since they consider a sce¬ 
nario outside the regime of applicability of the Carnot 
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FIG. 6: Quantum Otto Engine 



Schematic diagram of a quantum Otto engine, depicting only a single pair of levels, which can represent either 
a qubit, or the spacing of a harmonic oscillator. The four stages are (I) expansion from Ei to E 2 . (II) contact 
with a hot thermal reservoir. (Ill) compression from E 2 to E\. (IV) contact with a cold thermal reservoir. 


limit (much in the same way that a regular car engine, 
consuming fuel, does not violate the second law, since it 
is also outside the regime of applicability). Conversely, it 
is interesting that the net effect of squeezing appears to 
be as if the hot reservoir has been heated to a tempera¬ 
ture /3 h {r ), and that the performance of the machines is 
bounded exactly by the Carnot limit with respect to this 
new temperature. 


C. Quantum thermodynamic signatures 

One way to differentiate between a system which is 
genuinely using quantum effects and one which is only 
using the formal structure of quantum mechanics (the 
discreteness of energy levels, for example) is to devise 
signatures, or witnesses, for quantum behaviour. This is 
similar to what is done in entanglement theory, or in Bell 
nonlocality, where one finds witnesses which certify that 
entanglement was present, since no separable quantum 
state could pass a certain test. An interesting question is 
whether one can find analogous witnesses in a quantum 
thermodynamics setting. This is what was proposed in 
Ref. [234] in the form of Quantum thermo signatures. 

In more detail, the main idea of Ref. [234] is to find a 
threshold on the power of a thermal machine which would 
be impossible to achieve for a machine which is ‘classi¬ 
cal’. The authors take as the minimal set of requirements 
for a machine to be considered classical (i) that it’s oper¬ 
ation can be fully described using population dynamics 


(i.e. as a rate equation among the populations in the 
energy eigenbasis); (ii) that the energy level structure 
and coupling strengths are unaltered compared to quan¬ 
tum model under comparison; (iii) that no new sources of 
heat or work are introduced. A way to satisfy the above 
three constraints is to add pure de-phasing noise in the 
energy eigenbasis on top of the dissipative dynamics of 
the quantum model (arising from the interaction with 
the thermal reservoirs). One can then compare models 
with and without de-phasing noise, and ask whether the 
additional noise places an upper bound on the power of 
the machine. 

For simplicity in presentation, in what follows we will 
focus here on the results obtained for the four-stage qubit 
Otto heat engine, similar to the one described in the pre¬ 
vious subsection (except now with a qubit in place of a 
harmonic oscillator). We note that the authors show that 
the same results hold for a two-stage engine [235] and for 
continuous time engines [25], as well as for refrigerators 
and heat pumps. As an aside, the reason why the re¬ 
sult holds for all three models is because Ref. [234] also 
proves that in the regime of weak-coupling to the bath, 
and weak driving, all three types of engine can be shown 
to be formally equivalent, producing the same transient 
and steady state behaviour at the level of individual cy¬ 
cles. 

It is shown that a state independent bound can be 
placed on the power of a classical machine which is pro¬ 
portional to the duration of a single cycle of the engine 
T cyc , as long as the so-called ‘engine-action’ s is small, 
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where the engine action is the product of the duration 
r and energy scale (as measured by the operator norm 
of each term appearing in the Master equation). They 
demonstrate that there is a regime where a quantum en¬ 
gine (i.e. one without additional dephasing) can provably 
outperform the corresponding classical machine, with 
powers an order of magnitude larger in the former case. 


D. Stationary entanglement 

Entanglement is understood to be a fragile property 
of quantum states, that is one typically expects that 
noise will destroy the entanglement in a quantum state. 
Much effort has been invested in investigating and devis¬ 
ing ways in which one can counter the effects of noise, 
and maintain entanglement in a system, such as quan¬ 
tum error correction, dynamical decoupling, decoherence 
free subspaces, to name but a few. 

In the first subsection we saw that the non-equilibrium 
steady state of autonomous quantum thermal machines 
can be entangled. If one thus focuses not on their ther¬ 
modynamic functioning, but rather on their entangle¬ 
ment functioning, we see that whenever a thermal ma¬ 
chine reaches a steady state which is entangled, this con¬ 
stitutes a way of generating stead state entanglement, 
merely through dissipative interactions with a number of 
thermal environments at differing temperatures. 

Furthermore, if the interest is only in steady state en¬ 
tanglement generation, then it is not even necessary that 
the machine perform any standard thermodynamic task, 
and can in fact simply be a bridge between two reser¬ 
voirs, allowing the steady flow of heat from hot to cold 
such that the stationary state of the bridge is necessar¬ 
ily entangled. This is precisely the situation which was 
first considered in Ref. [236], where the minimal system 
of two qubits interacting with two baths at temperatures 
/?h and /3c was considered in the weak coupling (Marko¬ 
vian) regime. Numerous variants were then discussed: 
in Refs. [237-241] different aspects of the dynamical ap¬ 
proach to the steady state were analysed (assuming non- 
Markovian dynamics, the rotating wave approximation, 
etc); in Refs. [242, 243] a 3 qubit bridge was considered; 
in Ref. [244] the stationary discord was also studied; in 
Ref. [245, 246] geometric and dielectric properties of the 
environment were considered, and in Ref. [247] supercon¬ 
ducting flux qubits and semiconductor double quantum 
dot implementations were explored. 

Focusing on the simplest possible example, that of the 
two qubit bridge, the take home message of this line of 
investigation is that this is a viable means to generating 
stationary entanglement. In particular the implementa¬ 
tions considered in Ref. [247] suggest that in experimen¬ 
tally accessible situations steady state entanglement can 
indeed be maintained at a level which might be usable to 
then later distill. 


E. Outlook 

We have seen in this section a range of results con¬ 
cerning quantum thermal machines, focusing primarily 
on the quantum correlations and entanglement present 
in the machine, as well as other signatures of quantum¬ 
ness. Although we have focused on the progress that has 
been achieved so far, there are a number of directions 
which should be explored in further work to more fully 
understand the role of quantum information for quantum 
thermal machines. 

First of all, the main playground of study in this sec¬ 
tion has been the weak-coupling regime, where the ma¬ 
chine is in weak thermal contact with the thermal reser¬ 
voirs. It is important and interesting to ask what hap¬ 
pens outside of this regime, when the thermal baths are 
strongly coupled to the machine. On the one hand, in¬ 
tuition suggests that stronger coupling corresponds to 
more noise, which will be detrimental to fragile quantum 
correlations. On the other hand, stronger driving might 
lead to more pronounced effects. As such, the interplay 
between noise and driving needs to be better understood. 

Second, we have seen that quantum signatures, ei¬ 
ther in terms of entanglement or coherence, can be con¬ 
structed, which show that there is more to quantum ther¬ 
mal machines than just the discreteness of the energy 
levels. Here, it would be advantageous to have more ex¬ 
amples of quantum signatures, applicable in as wide a 
range of scenarios as possible. An experimental demon¬ 
stration of a quantum signature would also be a great 
development concerning the implementation of thermal 
machines. 

Finally, thinking of cooling as a form of error correc¬ 
tion, it is interesting to know if ideas from quantum ther¬ 
mal machines can be incorporated directly into quantum 
technologies as a way to fight de coherence. This would 
be as an alternative to standard quantum error correc¬ 
tion ideas, and an understanding of how they fit alongside 
each other could be beneficial from both perspectives. 

VII. FINAL REMARKS 

Ideas coming from quantum information theory have 
helped us understand questions, both fundamental and 
applied, about the thermodynamic behaviour of systems 
operating at and below the verge at which quantum ef¬ 
fects begin to proliferate. In this review we have given 
an overview of these insights. We have seen that they 
have been both in the form of technical contributions, 
for example with new mathematical tools for old prob¬ 
lems, such as the equilibration problem, and also in the 
form of conceptual contributions, like the resource theory 
approach to quantum thermodynamics. 

Although quantum information is only one of the many 
fields currently contributing to quantum thermodynam¬ 
ics, we expect its role to become more important as the 
field grows and matures. Indeed, we believe that plac- 
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ing information as a central concept, just as Maxwell did 
when his demon was born, will lead to a deeper under¬ 
standing of many active areas of physics research beyond 
quantum thermodynamics. 
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